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ABSTRACT 

We present a carefully vetted equatorial (± 30° Decl.) sample of all known single (within 
4") mid M-dwarfs (M2.5V-M8.0V) extending out to 10 pc; their proximity and low masses make 
them ideal targets for planet searches. For this sample of 58 stars, we provide Vj, Rkc, Ikc 
photometry, new low dispersion optical (6000 —9000 A) spectra from which uniform spectral types 
are determined, multi-epoch Ha equivalent widths, and gravity sensitive Nal indices. For 12 of 
these 58 stars, strict limits are placed on the presence of stellar and sub-stellar companions, based 
on a pioneering program described here that utilizes precise infrared radial velocities and optical 
astrometric measurements in an effort to search for Jupiter-mass, brown dwarf and stellar-mass 
companions. Our infrared radial velocity precision using CSHELL at NASA’s IRTF is ~90 m 
s“^ over timescales from 13 days to 5 years. With our spectroscopic results the mean companion 
masses that we rule out of existence are 1.5 Mjjjp or greater in 10 day orbital periods and 
7 Mji/p or greater in 100 day orbital periods. We use these spectra to determine rotational 
velocities and absolute radial velocities of these twelve stars. Our mean astrometric precision 
using RECONS [^(Research Consortium on Nearby Stars) data from 0.9-m telescope at Cerro 
Tololo Inter-American Observatory is ~3 milli-arcseconds over baselines ranging from 9 to 13 
years. With our astrometric results the mean companion masses that we rule out of existence 
are greater than 11.5 Mjjjp with an orbital period of 4 years and greater than 7.5 Mjjjp with 
an orbital period of 8 years. Although we do not detect companions around our sub-sample of 12 
stars, we demonstrate that our two techniques probe a regime that is commonly missed in other 
companion searches of late type stars. 


Subject headings: low mass stars, companions, planets 


1. Introduction 

Given human-kind’s search for life in the uni¬ 
verse, there is great motivation to find Earth-size 
and Earth-mass planets in the habitable zones 
of stars. Recent studies have determined that 
Earth size plane ts are c ommon around M-dwarfs. 
Morton fc SwiftI (|2013ll estimate an occurrence 
rate of 1.5 planets per M-dwarf with periods less 
than 90 days and radii larger than 0.5Kearth, 
using the list of 4000 stars with temperatures be¬ 


low 4000 K assembled in Batalha et al. ( 2013h 
that were observed with the Kepler mission 

’ I 2 OIOII . This 


( Borucki et al. 2010t Koch et ^ 


estimate is consistent with but slightly higher 
than previous studies, which measure occurance 
rates of approximately one planet per M-dwarf 


(Yolidin 20 111 Mann et al. 

2012: Swift et a.l. 20131: 

Dressine & CharbonnearJ 

2OI3II. Given the ap- 


parent abundance of Earth-size planets orbiting 
M-dwarfs, which dominate the stellar population 
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(75%: Henry et al.ll2006 ). Dressing fc Charbonneaul 
( 20131 ) predict the nearest non-transiting planet 


in the habitable zone orbiting an M-dwarf is 
within 5 pc of th Sun, with 95% confidence. 
However, how suitable these nearby planets 
within the classical ly defined habitable zone (e.g. 
Kasting et ah ^ 1993h may be for life is still under 


debate (e.g. iTarter et al. 2007 : Barnes et al ] l 20 nl 


ICnedel et alT 20l4l . Nevertheless, given their 

ubiquity and proximity, M-dwarfs are ideal stars 
to search for Earth-size and Earth-mass planets 
in stellar habitable zones. 

M-dwarfs have been favorite targets of precision 
searches for low mass planets, because a plane¬ 
tary companion will induce a greater reflex mo¬ 
tion on a low mass star than a Sun-like star, 
making it easier to detect. However, not all M- 
dwarfs are equally suitable targets for the precision 
measurements needed to find Earth-mass compan¬ 
ions. Some M-dwarfs have close stellar or substel- 
lar companions that may inhibit the detection of 
Earth-mass planets. The dynamically disruptive 
effects of these companions could also preclude the 
existence of Earth-mass planets; the lack of short 
period giant planets in mult iple planet systems 
corroborates this hypothesis dbatham et ah 2011 


Steffen et al.ll2012l : lBatalha et al.ll2013l) . Some M- 


stars also exhibit high chromospheric activity and 
large rotational velocitie s, which can hinder the 
achievable RV precision (iMohant v fc Basri _ 2003; 


Reiners fc BasrilEoOSH Jenkins et al. 2009t Reinerd 


2 OI 3 II : this is especially problematic for mid to late 
(M3 and cooler) M-dwarfs. Given the above con¬ 
siderations, we argue that the best stars to target 
for Earth-mass planet searches are likely the low¬ 
est mass stars (mid to late M-dwarfs) that do not 
have a disruptive companion, and are both inac¬ 
tive and slowly rotating. 

Yet, the statistics characterizing companions 
around mid to late M-dwarfs are still incomplete. 
Preliminary surveys show that Jupiter-mass com¬ 
panions are rare around M-dwarfs. Using ra¬ 
dial velocity ( RV) measurem e nts a nd high con¬ 
trast imaging, iMontet et al.l ( 20141 1 found that 
6.5+3.0% of M-dwarfs (M0-M5.5) host a gi¬ 
ant planet (1-13 Mj( 7 p) with a semi-major axis 
smaller than 20 AU, but this sample only included 
18 M-dwarfs of M4 or later in their survey of 111 
M-dw arfs. Another larg e M-dwarf survey of 102 
stars (jBonfils et al.1120131 1 only included M-dwarfs 


with V<14 and finds the giant planet (msini = 
100-1000 Mearth) frequency to be 2% for or¬ 
bital periods between 10 and 100 days, and the 
super-Earth (msini = 1-10 Mearth) frequency 
with orbital periods between 10 and 100 days to 
be significantly higher at 52%. Many of the M- 
dwarfs not surveyed are faint and chromospheri- 
cally active, which limits the achievable RV preci¬ 
sion and chances to find Earth-mass planets. As a 
result, these two surveys, which are representative 
of other spectroscopic M-dwarf surveys, include 
very few mid to late M-dwarfs and can report 
only preliminary statistics for planet occurrence 
around such stars. 

In order to obtain an unbiased assessment of the 
companion fraction of the nearest mid M-dwarfs, 
we construct a volume-limited sample out to 10 pc. 
Based on new, uniform optical spectra, we present 
this sample of 58 stars in §2. We report spectral 
types, Ha equivalent widths, and Na I indices for 
these stars in §3. We list our Vj, Rkc, IkS 
(hereafter without subscripts) photometry in §4. 
In §5, we focus on a sub-sample of 12 stars, for 
which we obtain high dispersion infrared spectro¬ 
scopic data described in §6 and astrometric data in 
§7. Results for the remaining stars in our volume- 
limited sample will be presented in a subsequent 
paper. We describe our Monte Carlo technique 
and rule out the existence of massive gas-giant 
companions, brown dwarfs and stellar companions 
in §8, and we conclude with a brief summary in §9. 


2. Sample Selection 


2.1. An Equatorial Sample of Nearby Mid 
M-Dwarfs 

Beginning with a parent volume-limited sam¬ 
ple of stars extending out to 10 pc (unpublished 
list; iHenrv et al.l 120061 1. we assembled a sample 
of mid M-dwarfs for detailed study outlined in 
Table [TJ These stars are within the declination 
range of + 30°and are thus accessible to the ma¬ 
jority of observing facilities in both the north¬ 
ern and southern hemispheres. To be inclusive, 
we define mid M-stars using three independent 
measures, including optical spectral type, (V-K) 
color and absolute magnitude. Starting with the 


^Subscripts: J indicates Johnson and KC indicates Kron- 
Cousins. 
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complete 10 pc sample, we include stars meeting 
one of the following criteria: spectral types be¬ 
tween M3.5V and M8.0V classified by the RE¬ 
CONS system (described in §3), V-K =5.0—9.0 
mag or My = 12.0—19.0 mag. There are 69 systems 
that meet at least one of these criterisd. Of those 
systems, 13 are known binaries, five are known 
triples, one is a quintuple system (GJ 644ABCD- 
GJ 6 431 and one is a kn own multi-planet host (GJ 
876; iRivera et al. The 21 close-separation 

(<4") mid M-dwarf binarie^ within the distance 
and declination range of this sample are listed in 
Table [2] Also listed are 2MASS coordinates, par¬ 
allaxes, spectral types, absolute V magnitudes, V, 
A, / apparent magnitudes, near infrared photom¬ 
etry from 2MASS (J, H, Kg apparent magnitudes) 
and the configuration of the system. 


A sub-sample of mid M-dwarfs is constructed 
that excludes these close binaries (<4") and the 
planetary host GJ 876. This sub-sample includes 7 
mid M-dwarfs that are wide companions to higher 
mass stars (GJ 105B, GJ 166C, GJ 283B, GJ 
644C, GJ 752B, GJ 896B, GJ 1230B). At the start 
of this program, 60 mid M-dwarfs met the sample 
requirements. Subsequently, GJ 867B has been 
determined to be a single-lined sp ectroscopic bi¬ 


nary with a period of 1.795 days (iDavison et al 


I 2 OI 4 I I . Likewise, LHS 1 610 has been claim ed to be 
a spectroscopic binary ( Bonfils et al.ll2013li . but its 
orbital properties are unknown. Table [I] lists the 
astrometric, photometric, spectroscopic and phys¬ 
ical properties of the remaining 58 stars, which 
we refer to as effectively single equatorial mid M- 
dwarfs. Mass estimates in T able [D are based on 
mass lumii iositv relatio n s of iHenrv fc McGarthv 


(|1993D and iHenrv et al.l (|l999l) : typical errors us¬ 


ing these relations are close to 20%. Because the 
long term goal of this program is to conduct a more 
comprehensive 3D search for companion to these 
stars, we refer to this sample as GAESAR, which 
stands for a Gompanion Assessment of Equatorial 
Stars with Astrometry and Radial velocity. 


Of the 58 stars identified above as possible tar¬ 
gets for precision planet searches among equato- 


®We include GJ 628 with V-K=4.99 in our sample, as its 
V-K value may be greater than 5.0, given our photometric 
uncertainties. 

^This list includes four mid M-dwarf binaries that are wide 
companions to more massive stars (LP 771-096BC, GJ 
569BC, GJ 695 BC, GJ 867BD). 


rial mid M-dwarfs, only 33 of the stars have been 
included in past spectroscopic searches for plane¬ 
tary companions dBarnes et al 1126121: iRodler et al 


2012 : Tanner et all 2012 : Bonfils et a,l.l2013r Montetetal 


20141) . In addition, four of the 58 stars do not have 
known projected rotational velocity (wsini) values, 
and thus could be rapidly rotating (see Table [T|) . 


3. Optical Spectroscopic Measurements 
3.1. Observations 

We obtained optical (6000 — 9000 A) spec¬ 
tra of all 58 stars in the CAESAR sample be¬ 
tween 2003—2006 and 2009—2011 using the Cerro 
Tololo Inter-American Observatory (CTIO) 1.5-m 
Richey-Chretien Spectrograph (RCSpec) with the 
Loral 1200 x 800 CCD camera, as part of the 
broader RECONS spectral-typing program (e.g. 
Henry et al. 2004t Jao et all 2008). The spectra 
were obtained with the #32 grating in first order 
with a 15.1° tilt, which yields a spectral resolu¬ 
tion of 8.6 A; the spectra were acquired through 
the OG570 order blocking filter. For consistency 
checks and to mitigate the effects of cosmic rays, 
two spectra of each target were taken consecu¬ 
tively. In addition, the majority of stars have 
spectroscopic observations on multiple epochs. 
To assist with spectral classification, at least one 
flux standard was observed each night, and an 
ensemble co l lectio n of spectral standards from 
Henry et al.l (120021) were observed. Some of the 


stars in the GAESAR sample are used as stan¬ 
dards, namely GJ 283B, GJ 644G, GJ 752B, GJ 
1065, GJ nil, GJ 1154, GJ 1207, LHS 292, LHS 
2090 and LHS 3799. 

The data were reduced with standard IRAF 
techniques; bias subtraction and dome and/or sky 
flat holding were performed on the data using cal¬ 
ibration frames taken at the beginning of each 
night. Fringing was effectively removed from the 
data using a combination of dome and sky flats. 
One flux standard per night was used for abso¬ 
lute flux calibration. Spectra were wavelength 
calibrated using consecutively recorded HeAr arc 
spectra. Further detai ls regarding r educt ion and 
extraction are given in IHenrv et al. ( 2004 ). 


3 


















































3.2. Results 

3.2.1. Spectral Types 


To assign a spectral type, the wavelength cal¬ 
ibrated spectra are resampled via interpolation 
onto a fixed 1A grid. The Ha and tell uric features 


based on the sky transmission map of lHinkle et al 


( 2 OO 3 I I are then given a spectral weight of zero to 
essentially remove these features from our analy¬ 
sis. The spectra are normalized to a value of 1 at 
7500 A. Then, the spectra are then compared to 
the library of observed standards (see Figure [U, 
and the adopted spectral type is that of the stan¬ 
dard that yields the lowest standard deviation 
of the target spectrum divided by the standard 
spectrum over the spectral range (6000 — 9000 A). 
The determined spectral types range from M2.5 to 
M8.0, and are listed in Table [T] The uncertainty 
in all cases is ±0.5 spectral sub-classes, based on 
consistency between multiple epochs. Additional 
details o f the spectral type d eterminations are pro¬ 


vided in Riedel et al. ( 2014l l. 


3.2.2. Ha Emission 

To assess the amount of chromospheric activ¬ 
ity that these stars exhibit, we measure equiv¬ 
alent widths of Ha 6563 A. To do this, we first 
subtract the linear continuum from the region be¬ 
tween 6550 A and 6575 A and then assume the 
strongest maximum (or minimum) within this re¬ 
gion is the Ha line. After determining its loca¬ 
tion, we integrate both the continuum and spec¬ 
trum over 22 A, with the central wavelength cor¬ 
responding to the Ha line to determine the equiv¬ 
alent width. Upon visual inspection of the spec¬ 
tra, we could not reliably determine emission or 
absorption lines smaller than 0.5 A. Therefore, we 
conservatively use this number as the error value 
on the Ha measurements, which are given in Ta¬ 
ble 131 we adopt the standard convention of de¬ 
noting emission with a negative sign. If any star 
exhibits Ha equivalent widths less than the the 
typical noise level of 0.5 A in at least one epoch, 
we categorized it as an emission-line star. This is 
denoted with an ‘e’ next to the star’s spectral type 
in Table [T] 

We classify 35/58 of the stars from the CAE¬ 
SAR sample as emission-line stars. We also note 
that the fraction of emission-line stars increases 
with later spectral type. 37% of the 35 stars with 


spectral types M4 or earlier are emission-line stars, 
whereas 100% of the 9 stars with spectral type 
M6 or later are emission-line stars. Our fraction 
of earl y M-dwarfs with emis sion is comparable to 
that of iReiners et al. ( 2012 1 who find 41% of the 
115 stars M3 to M4.5 to b e activ e emission-line 
stars. Likewise, Gizis et al. ( 2000l l reports an in¬ 
crease in emission-line stars with lower mass and 
finds 100% of the M7 stars are emission-line stars. 
He also notes that this trend breaks down past M7 
and that the fraction of later type stars that show 
chromospheric activity is significantly less. 


3.2.3. Surface Gravity Indices 

To assess the surface gravity of these stars, 
as a possible tracer of their evolutionary states, 
the gravity sen sitive Nal doublet (defined in 
Lvo et all 2004l l is measured; this feature is well 
known to be weak in giants and strong in dwarfs 


__ _ ng 

(^e.g. Allers et al. 2007; Schlieder et ahllioi^ . For 
the Nal doublet index, we use the wavelength 
region between 8148A and 8200A. We conserva¬ 
tively estimate our index error to be 0.05; this 
value is the average difference value of our mea¬ 
sured Na I doublet indices for all of our stars with 
more than one epoch. These values are given in 
Table [3l While these indices derived from the low- 
resolution spectra are a good indicator of the evo¬ 
lutionary state of the star, high-resolution spectra 
are needed to ascertain a more quantitative mea¬ 
sure of the surface gravity. Also, we do caution 
that these lines can be affected by metallicity and 
stellar activity; metal poor stars and chromospher- 
ically active stars will syste matically have lowe r 
equivalent widths (EW) fe.g. iHawlev et al.] 1999[l . 

All of the stars in the ensemble sample ex¬ 
hibit average Nal doublet indices in a typical 
range for main sequen ce stars (>1.1; ILvo et ^ 


2004 lAllers et al.l [20071) ; we therefore classify all 
58 CAESAR stars as on the main sequence. This 
is denoted by adding a ‘V’ to the spectral type 
listed in Tabled] 


4. Optical Photometry 

Prior to our measurements, four of the 58 stars 
from CAESAR sample did not have complete sets 
of V, R and I photometry. The remaining stars 
had photometric measurements presented in eight 
different publications. Therefore, to establish a 
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uniform, homogenous, set of photometric measure¬ 
ments for the ensemble sample of 58 stars, we ob¬ 
tained optical photometric observations using the 
0.9-m telescope at CTIO. For all of our photom¬ 
etry frames, we use the center 1024x1024 pixels 
on the Tektronix 2048x2048 CCD. The CCD chip 
has a plate scale of 0'.'401 pixel“^, which gives 
a field of view (FOV) of 6.8 by 6.8 arcminutes 
( Jao et aP 2003l l. All frames were collected at an 
airmass less than 2 and with the target star having 
a signal-to-noise ratio (S NR) > 100. We use 10 or 
more standa rd stars f rom Landolt ( 1992 ), iLandoltl 
(I 2 OO 7 II . and ICra.ha.ml (|l982h to create extinction 
curves each night, and transformation equations 
to obtain V, R and I photometry for all our target 
stars and reference stars used for astrometric mea¬ 
surements described in §7.2. During the course of 
observations, we used two dif ferent V filters for 
photometry. I.Tao et all (201 ll ) demonstrate that 
both V Hlters give effectively identical V band 
photometry for standard stars; therefore it is suit¬ 
able for us to combine photometry from the two 
filters. For additional details on the photometric 


reduct ion and its associated err ors, see I.Tao et al 
( 2005 1 and Winters et al.l ( 201ll ). 

The V, i?, / photometry and the number of 
nights on which observations were made are re¬ 
ported in Table [TJ Errors at V, R and / are 
0.02—0.03 mag. A comparison of five stars (GJ 
300, GJ 406 GJ 555, GJ 628 and GJ 729) with 
Bessel ( 1990l) indicates that the two datasets are 


consistent to 0.04 mag. Using our photometric 
measurements and the parallax data given in Ta¬ 
ble [H the absolute magnitude errors range from 
0.03 to 0.08 mag. 


5. A Companion Search of a 12 Star Sub¬ 
set 

We present results of a 3D companion search 
on a sub-sample of 12 CAESAR stars, including 
G 99-49, GJ 300, GJ 406, GJ 555, GJ 628, GJ 
729, GJ 1002, GJ 1065, GJ 1224, GJ 1286, LHS 
1723, and LHS 3799, which are the most data-rich 
in our sample. The remaining stars in the en¬ 
semble sample will be presented in a subsequent 
paper. These twelve stars have astrometry base¬ 
lines, ranging from 9 to over 13 years, and have at 
least 5 infrared radial velocity measurements span¬ 
ning from almost 2 weeks to 5 years. The spectral 


classes of these stars range from M3.0 to M5.0. 


6. High Dispersion Infrared Spectroscopy 
6.1. Observations 

All infrared spectroscopic observations were 


Ail mirarea spectrosco pic ooservations were 
obtained using CS HELL (ITokunag-a et al. 1990l: 


Greene et al.lll99,3ri located on the 3-m telescope 


at NASA’s Infrared Telescope Facility (IRTF). 
CSHELL is a long-slit echelle spectrograph that 
uses a circular variable filter to isolate a single 
order onto a 256x256 InSb detector. Spectra are 
centered at 2.298 microns (vacuum) and cover 
approximately a 50 A window. Telluric methane 
absorption features from the Earth’s atmosphere 
that are superimposed on the photospheric ^^CO 
R branch lines at 2.3 microns a re used as an abso¬ 
lute wavelength re ference (e.g. Blake et ^ 2010l: 
Bailey et al. 20121) . The design resolving power 


is 100,000 per pixel. We use CSHELL in the 
high resolution mode (0.5" slit=2.5 pixels), which 
yields a predicted spectral resolving power of R 
~ 40,000. The measured resolving power deter¬ 
mined from the model fits to telluric absorption 
features (described later) is ~57,000. This number 
is significantly higher than the predicted resolv¬ 
ing power for CSHELL using the 0.5" slit. We 
note this discrepancy because CSHELL has an 
adjustable slit; it may be that the slit is s maller 
tha n the designed 0.5 ". Crockett et al.l ( 2011 1 
and IPrato et al. ( 20081) also report a higher than 
predicted spectral resolving power (R~46,000) 
for CSHELL in high resolution mode. Also, we 
note a similar effect of determining a higher spec- 
tral resolving pow er for Keck is determined in 
iBailev et al. ( 2012 1. and may be a feature of the 
analysis code. 

Each night, two spectra of a given star were 
obtained in succession at two different positions 
along the slit, separated by 10". Hereafter, we re¬ 
fer to these two positions as nod A and nod B. 
Observations were obtained between 2008 Novem¬ 
ber and 2014 January. 

Our exposure times ranged from 180 to 1200 
seconds per nod position, and were set to yield 
SNRs of 125 per pixel (or optimally, a combined 
SNR of 175-I-) for most of our targets. For the 
faintest three stars {Kg > 8.0) in our sample of 12 
stars (GJ 1065, GJ 1224, GJ 1286), a SNR of 125 
per pixel was not achieved as the maximum expo- 
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sure time is set to 1200 seconds to limit cosmic ray 
events and the dark current. 

At the beginning of each night, we obtained 
a minimum of 30 flat and 30 dark images each 
with an integration time of ten seconds. Also, 
on November 14, 2009, we collected an additional 
30 flat and 30 dark images with an integration 
time of 20 seconds, which were used in creating a 
bad pixel mask described in §6.2. Most nights we 
also observed bright stars of spectral type early 
A, as these stars exhibit no intrinsic absorption 
or emission lines in this wavelength region and 
therefore can be used to identify telluric features. 
These telluric standards are used to characterize 
the instrumental profile and wavelength solution. 
When first collecting the data, we did not real¬ 
ize how sensitive our final RV measurements were 
to the initial solutions for the instrumental pro¬ 
file and wavelength solution obtained from the tel¬ 
luric standards. After reducing part of the data, 
we determined that A star observations obtained 
nightly yield the best precision. In a few cases, we 
only obtained a few A star observations per run 
leaving us with four nights that contain no A star 
observations. The four nights the telluric stan¬ 
dards were not observed are marked in Table U by 
an asterisk next to the date. 

We aimed to observe each target at least four 
nights within a single observing run (~ 1—3 weeks) 
in order to search for companions with orbital pe¬ 
riods of less than a week. Because of inclement 
weather we were not able to achieve this ca¬ 
dence for all targets. On subsequent runs, we re¬ 
observed our targets at least once to search for 
companions with longer orbital periods, except 
for GJ 1065. For the sample of 12 stars studied 
here, we obtained between 5—12 RV epochs for 
each star, spanning a temporal baseline between 
13 days and 1884 days (see Table S]). 

6.2. Image Reduction and Spectral Ex¬ 
traction 

We subtracted each nodded pair of images from 
one another to remove sky emission, dark current 
and detector bias assuming that changes in the 
detector or spectrograph properties were negligi¬ 
ble over the timescale when the nodded pair of 
images were obtained. After completing the nod- 
subtraction, we corrected each image for flat field¬ 
ing. Corrections for flat fielding were performed by 


generating a nightly master flat field image from 
all flat field images obtained on a particular night. 
The master flat held images were created by hrst 
subtracting the median dark image of the same 
exposure time from each of the hat held images. 
Then, each hat held image was normalized to the 
central 15% of the array, which was the bright¬ 
est section of the array and the least affected by 
deviant pixels. After normalizing the image, all 
images were median combined. 

We then applied a bad pixel mask to our spectra 
to remove dead and hot pixels from the data. To 
identify dead pixels, we located any pixel hve times 
below the standard deviation of the median pixel 
value of the master hat held array. To locate hot 
pixels, we subtracted two times the count value 
of the 10 second exposure master hat held image 
from the 20 second exposure master hat held im¬ 
age, and then normalized this number to the 20 
second exposure master hat held. Because these 
pixels should increase linearly with time and there¬ 
fore have the same values, we identihed any pixel 
with values greater than three times the standard 
deviation of this median pixel value of the differ¬ 
ence image to be a hot pixel. All deviant pixels 
identihed from the bad pixel mask are assigned 
interpolated values using the neighboring pixels. 


We optimally extract ed each spect rum follow¬ 


ing the procedures in Horne ( 1986ll as imple¬ 


mented for nod-subtracted spectra in iBailev et al 


(|2012i) . The code used to analyze the data in this 
work is a modihed version of that described in 
Bailey et al.l ( 2012l l. tuned to work for CSHELL 


data. The advantage of optimally extracting spec¬ 
tra over the standard extraction is that the opti¬ 
mal extraction minimizes the noisy contribution of 
the prohle wings and eliminates and/or mitigates 
noise features within the spectral prohle caused by 
cosmic rays and deviant pixels not excluded with 
the bad pixel mask. 


To obtain an optimally extracted spectrum, we 
summed the pixels from the nod subtracted images 
over the cross-dispersion to give the standard spec¬ 
trum. We then htted a second order polynomial 
to map the curvature of the order on the detector. 
Next, we htted the spectral prohle of the standard 
spectrum with a Gaussian to model our spatial 
prohles at each pixel step (column) along the or¬ 
der parallel to the dispersion direction of the nod- 
subtracted spectral image. From this, the variance 
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of the profile was determined. Then, we summed 
the pixels weighted by the variance image of the 
spectrum’s spatial profile to create our two dimen¬ 
sional optimally extracted spectrum. For low SNR 
data, we implemented a clipping routine that in¬ 
terpolates over a pixel that is more than 5cr above 
or below the running average of the five pixels next 
to the pixel in question to remove any remaining 
deviant features. 


To obtain an estimate of the SNR for each 
spectrum, we use a simplified ver s ion of the CCD 
equation from Mortara fc Fowler ( 198l[ l modified 
to account for the noise introduced by subtract¬ 
ing pairs of images. When performing a nod 
pair subtraction we remove the sky background, 
dark current and bias simultaneously. Therefore, 
we cannot distinguish between these values and 
refer to them collectively as the uncertainty in 
background. Executing a pair subtraction means 
that we have to deal with this uncertainty in the 
background twice and the read noise associated 
with each of those background estimates. In most 
cases, the background of the first image (Image 
A) should equal (within the uncertainty) the back¬ 
ground of the second image (Image B). Therefore, 
we simply double the noise contribution from the 
background and the read noise. The equation to 
calculate the SNR per pixel is as follows: 


SNR= , 

^ySe + 2■n{Be+Rl) 

where Se is the total number of counts per inte¬ 
grated column of a spectrum extracted from a nod- 
subtracted image in electrons, n is the number of 
pixels in the spatial direction that are integrated 
over during the extraction. Be is the integrated 
background counts of the corresponding sky image 
before image subtraction is performed, in electrons 
per pixel, and Re^ i s the read no i se set to be 30 
electrons/pixel from iGreene et al.l (|l993r) . 

Following the above description, we determine 
the SNR for each integrated pixel of the spectrum 
and set the final SNR for the spectrum to be the 
mean of these values. The SNRs are then added in 
quadrature for the nod A and nod B measurements 
to give a combined SNR value. 

Because of occasional poor weather conditions 
leading to low SNR, not all observations are suit¬ 
able for precision RV analysis. We require the 


SNR of the individual spectrum in the nod pair 
to be greater than 50 and the reduced estimate 
of our modeling prescription (§6.3) to be below 

3.5. 


6.3. Method to Determine Spectral Prop¬ 
erties 


We fit each observation to high resolution spec¬ 
tral models that are convolved to the resolution 
of CSHELL. Each model spectrum is formed by 
combining a synthetic stellar spectrum and an 
empirical telluric spectrum. The synthetic stel- 
l ar spectra are create d from NextGen models 
( Hauschildt et al.lll999ll . The telluric model spec¬ 
tra are extracted from observations of the Sun 
from an ul tra high resolution KPNO/F TS telluric 
spectrum ( Livingston fc Wallacelll991 ). We adopt 
the stellar template closest in temperature to our 
star, using the assi gned spectral types and th e 
temperature scale of Kraus fc Hillenbrand (l2007ll . 
We fix the surface gravity log(g) to 4.8 dex (cgs) 
for all of our stars, whi ch is consistent with mea¬ 
surements assembled in Mentuch et al. ( 2008^ and 
Hillenbrand fc White! ( 2004 ) for field M-dwarfs. 


The model spectrum consists of 19 free parame¬ 
ters to fit. The linear limb darkening coefficient is 
set to 0.6 for all stars, which is appropriate for cool 
stars at infrared wavelengths ( Glaretll2000ri . Three 
of the parameters make up a quadratic polynomial 
that characterizes the wavelength solution. Nine 
of the parameters are Gaussians used to model the 
line spread function (LSF) of the spectrum; we 
assume that the LSF along the order is constant. 
The remaining six parameters are the depth of the 
telluric features, the depth of the stellar features, 
the projected rotational velocity (usini), the RV, 
a normalization constant, and a linear normaliza¬ 
tion term. 


We fit the empirical telluric spectrum to our 
rapidly rotating A star for each night we pro¬ 
cured observations of A stars. From this measure¬ 
ment, we estimate the wavelength solution and 
the instrumental profile. The instrumental pro¬ 
file is solved for by interpolating the input spec¬ 
trum onto a log-linear wavelength grid and con¬ 
volving it with a Gaussian kernel set to the spec¬ 
tral resolution of the instrument determined by 
fitting the telluric spectrum. We tested both sin¬ 
gle and multiple Gaussian functions to obtain the 
instrumental profile of GSHELL. We favor multi- 


7 










































pie Gaussians to fit the instrumental profile, as we 
had better agreement between the RV estimates 
from the AB nod pairs and smaller RV dispersions 
overall. The Gaussian kernel is composed of one 
central Gaussian and four satellite Gaussians on 
ea ch side followi n g clos ely the technique described 
in IValenti et al.l 


(I 995 I). We set the positions of 
the centers of the Gaussians and the widths of 
the satellite Gaussians so that the curves barely 
overlap. The amplitude of the central Gaussian is 
constrained by the normalization factor, while the 
width of the central Gaussian is allowed to vary. 
The amplitudes of the satellite Gaussians are al¬ 
lowed to vary. Optimization of these values is ac¬ 
complished by minimizing the varia nce weighted 


reduce d chi-squared as described in I Bade v et al 


( 2 OI 2 II to best reproduce the observed spectrum. 
In the cases when no A stars were observed on 
a night, we use the mean values determined on 
nights close to the night when A stars were ob¬ 
served. 

After fitting the telluric spectrum, the nine pa¬ 
rameters used to characterize the LSF are kept 
constant for all remaining fits. We use an itera¬ 
tive process where we fit the target spectrum to 
the combined synthetic stellar model and empiri¬ 
cal telluric model. On the first iteration, we fit the 
wavelength solution, the depth of the telluric spec¬ 
trum, the RV, the normalization constant, and 
the linear normalization term. With an improved 
guess on our second iteration, we allow the usim, 
the depth of the telluric model, the depth of spec¬ 
tral model and the two normalization constants to 
fluctuate. The vsini is determined following the 
description provided in Gray (20^. We adopt 
the average usini value from this iteration for all 
epochs as the nsini value for the star. Finally, 
we repeat the hrst iterative process allowing the 
wavelength solution, the depth of the telluric spec¬ 
trum, the RV, the normalization constant, and the 
linear normalization term to vary in order to de¬ 
termine the absolute RV of the star. Computa¬ 
tionally, the optimization of the model spectrum 
is completed using AMOEBA, which is a routine 
used for minimization of multiple variables using 
the d ownhill simplex method of iNelder fc Mead 
(1965). We note that AMOEBA is very sensi¬ 
tive to initial guesses and is given user specified 
ranges to restrict the answers to physically rea¬ 
sonable solutions. An example of an optimally 


extracted spectrum fit to our telluric and stellar 
models is shown in Figure [2] 

Rather than use the full 256 pixels along the 
order, the modeling analysis is restricted to pixels 
between 10 and 245, which corresponds to a small 
continuum area on the spectra. These boundaries 
are set to prevent strong absorption features from 
moving in and out of the analysis region on differ¬ 
ent epochs, because of different barycentric correc¬ 
tions. Partial features that are cut off by the edge 
of the chip can cause our RV value to change on 
the order of 100 m s“^. Using the restricted pixel 
range, our average precision improved by 27% for 
the 12 stars analyzed here. 


6.4. Spectroscopic Results 


The RV results of the spectroscopic modeling 
are given in Table |4] with objects listed alphabeti¬ 
cally. Multiple RV measurements on a single night 
are averaged to provide a single epoch value. All 
RV measurements are corrected to the Solar Sys¬ 
tem Barycenter using a correction prescription ac¬ 
curate to ~ 1 m s“^ (G. Basri; priv. communica¬ 
tion) . 


Under the assumption that the stars do not 
have companions, the observed dispersion is 
thought to be caused by a combination of theoret¬ 
ical photon noise error, intrinsic stellar error and 
instrumental error (aobs'^ = cFphoton'^ + cTsteiiar'^ 
+ cfinstr^)- The theoretical error for each spec- 
trum is calculated based on the prescription by 
Butler et al. ( 1996ll . For our 12 stars, the average 
theoretical photon noise error is 57 m s“^, with a 
standard deviation of 14 m s“^. The stellar error 
is assumed to be zero, as we are observing a field 
population of slowly r otating M-dwarfs. This is 
supported by results of Bonfils et al. ( 201,111 show¬ 
ing that the observed RV dispersions based on 
optical spectral for nine of the stars in this sub¬ 
sample are below 10 m s“^. We solve the equation 
above to determine the instrumental error {cfinstr^ 
= <^obs^ — O’photon^) for each star, in which the ob¬ 
served dispersion is the standard deviation of the 
nightly radial velocity measurements and the the¬ 
oretical photon error for each star is the average 
of the nightly theoretical photon errors. The aver¬ 
age instrumental error for this subset is 73 m s“^, 
with a standard deviation of 42 m s“^. We adopt 
this number as the instrumental error for all tar¬ 
gets in our sub-sample. In Table IH we report the 
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Fig. 1.— CTIO spectra of five CAESAR stars 
(black) normalized to 7500 A. The comparison 
spect ra (red) are from the library of standards 
from iHenrv et ah ( 2002 ). Small discrepancies at 
the longest and shortest wavelengths are due to 
minor errors in the flux calibrations. The vertical 
line represents the position of Ha (6563 A). 
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Fig. 2.— Spectral Modeling of GJ 300. Spectra 
are modeled by combining a telluric spectrum (top 
spectrum) with a synthetic stellar spectrum (2"^^ 
spectrum); the telluric spectrum provides an abso¬ 
lute wavelength reference. The CSHELL spectra 
of GJ 300 are shown (black) in comparison with 
the the best fit (red; 3’’'^ spectrum). The residuals 
of the fit are shown (bottom spectrum). 
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final error assigned to each measurement, which 
is calculated as the instrumental error added in 
quadrature with the theoretical photon noise er¬ 
ror. 


In Table [5l we summarize the infrared spec¬ 
troscopic results, including the absolute RV, the 
number of epochs, the time span of observations, 
the standard deviation of the RV measurements 
and the nsini value and its uncertainty. The abso¬ 
lute RV is the mean of the RV measurements from 
different nights. We note that systematic uncer¬ 
tainties in the adopted synthetic template (log(g), 
Te//), and the wavelength region used in the fit 
can cause RV shifts of ~ 100 m s“^. Therefore, 
we set the uncertainty of the absolute RV measure¬ 
ments to be 100 m s“^ for all 12 stars. All of our 
stars have previous absolute radial velocity mea¬ 
surements and those measurem ents are less than 
3 sigma from our m easurements ( Gizis et al.ll2002 : 


Nidever et al.ll2002D . 


The vsini value is the average of the nightly 
best fit usinf measurements. The error on the 
vsini value is calculated as the standard devia¬ 
tion of the best fit nightly vsini measurements. 
We do caution that the spectral resolving power 
of CSHELL is not high enough to fully resolve the 
lines of the slowest rotators. Line broadening be¬ 
comes measurable for vsini values in excess of 3 
km s“^, therefore we set this value as our nsini 
detection limit. This de tection threshold is i n line 
with those reported by Reiners et al.l (1201211 of 3 
km s“^ and Browning et al. I (12011)11 of 2.5 km s"! 
for similar resolution spectra (R=45,000-48,000). 
We detect rotational broadening above our detec¬ 
tion threshold for two stars, G 99-49 and GJ 729, 
out of the twelve. The prev ious r>sinf value for G 
99-49 of 7.4+0.8 km s“^ by Delfosse et al. ( 1998h 
is within 2 sigma of our measurement of 5.8+0.3 
km s“^. Likewise, the pre vious nsinz value for GJ 
729 of 4.0+0.3 km s“^ bv iBrowning et idl ( 20101 ) 
is within 0.7 sigma of our measurement of 3.8+0.6 
km s“^. 


Our observed dispersions for these 12 stars 
range from 47 m s“^ to 139 m s“^. Our average 
observed dispersion is 99 m s“^, with a standard 
deviation of 27 m s“^. Figure [3] shows the RV 
measurements for all epochs for each of the target 
stars. The median RV error for high SNR spec¬ 
tra of the 10 slowly rotating stars (?;sim < 3.0 km 
s“^) is 88 m s“^. This is dominated by the instru- 
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, which suggests a limiting 
for high SNR, slowly ro- 


mental error (73 m s“^ 
precision of ~90 m s 
tating mid to late M-dwarfs. We note that some 
of the instrumental uncertainty could be a conse¬ 
quence of our modeling prescription; a more so¬ 
phisticated approach may yield better results. We 
also note that a precision of 58 m s“^ has been 
reported for multiple epoch meas urements of the 


MO s tar GJ 281 using CSHELL (jCrockett et al 


20111). These precisions are nevertheless consid¬ 


erably better thy the design specs for C SHELL 
( Tokunaga et al.iri99fll : Greene et al.lll993h . espe¬ 
cially considering the small wavelength coverage, 
and are credited to the talents of the instrument 
team. 


7. Astrometry 
7.1. Observations 

All optical astrometric observations were made 
using the 0.9-m telescope at CTIO. The astrome¬ 
try program began as an NOAO Surveys Program 
in 1999 August and continued from 2003 Febru¬ 
ary as part of the SMARTS (Small and Moder¬ 
ate Aperture Research Telescope System) Consor¬ 
tium. Stars have been intermittently added to the 
observing list since 1999 and stars discussed here 
continue to be observed. The Cerro Tololo In- 
teramerican Observatory Parallax Investigation 
(CTIOPI) program was originally designed to 
measure accurate parallaxes of nearby stars. We 
are now using the same data and techniques to 
look for perturbations that remain in our astro- 
metric signal after solving for the parallactic mo¬ 
tion and the proper motion of our targets. The 
presence of a periodic perturbation in our data 
might signify that our star is orbiting around a 
common center of mass with an unseen compan¬ 
ion. To do this, we use the same instrumental 
setup as that used for the photometry frames. We 
observe each star in one of the V, R, or / filters. 
Stars are observed through the hlter that gives the 
strongest reference field, while not compromising 
the counts of the target star (the filter and num¬ 
ber of reference stars used are given in Table El). 
Strong reference fields that give the most precise 
parallax measurement include 5 to 12 reference 
stars that are bright (peak counts greater than 
1000 ), close on the chip to the target star, and 
in a configuration that surrounds the target. We 


require all stars used as reference stars to have a 
minimum of 100 counts. Our reduction routine 
accounts for plate scaling and rotation effects, but 
ignores higher-order ter ms (astigmatism , coma, 
chromatic aberration; see Ijao et aLll2005ll . 


Exposure times are set such that the target star 
or in some cases a very close bright reference star 
does not saturate. Our maximum exposure time 
is 600 seconds. We aim to obtain exposure times 
of at least 30 seconds for every star, although this 
is not always possible for our brightest stars. 

Frames are only collected under seeing condi¬ 
tions better than 2.4", determined by the FWHM 
of the stars in the field to be used in our reduc¬ 
tion. Also, the target star and its reference stars 
must have an ellipticity less than 20%, in order 
to determine the centroid of the stars with the 
most accuracy and to eliminate frames with pos¬ 
sible tracking/guiding errors. The guider is typi¬ 
cally used for any exposure times longer than 300 
seconds. 


Eight of the twelve stars discussed here are ob¬ 
served astrometrically with the V filter. During 
the course of observations, the first V filter was 
cracked and replaced by another V filter. The use 
of the second V filter from 2005 February to 2009 
July causes a few milli-arcseconds (mas) offset in 
astrometric residuals of known singles from other 
techniques. In 2009 July, we switched back to the 
original V filter, as the minor crack on the filter 
edge does not affect the data acquired on our cen¬ 
tral quarter region of the GGD. Using data from 
both hlters gives the same paral lax measurement 


but with slightly higher errors (jSubasavage et al 


20091 : Riedel et al. 2010l) and the average residual 
deviation for the stars is still less than 4 mas for 
our 12 stars (see §7.3). Therefore, we choose to 
use data obtained in both filters to maximize the 
time coverage. 


7.2. Astrometric Reductions 

We correct all centroids of reference stars and 
the target s tar for differen tial chromatic refrac¬ 
tion (DCR; Ijao et ah 20051) . We measure accu¬ 
rate positions using the SExtractor Ge ntroiding 
algorithm from Bertin fc Arnouts (Il996l) a nd use 
the Gaussht program ( Jeffervs et al.l 198^ to si¬ 
multaneously solve for the parallax relative to the 
reference stars and proper motion on all available 
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data (for more details see Jao et al. 20051) . If after 
running our Gaussfit program, reference stars are 
found to have proper motions greater than Ol'05/yr 
or a parallax greater than 5 mas based on photo¬ 
metric parallaxes, then those stars are rejected as 
reference stars. To obtain the absolute parallax 
of our target star, we must correct for the paral¬ 
lactic motion of the references stars as these stars 
are not infinitely far away. We use photometric 
parallaxes and accurate V, R and I photometry 
described in §4 to correct o ur relative parall ax to 
the absolute parallax value ( Jao et ablbnObl) . All 
the frames are used to fit the parallactic orbit to 
the star. However, only nights with 2 or more 
good images are used in order to calculate the as¬ 
trometric signal that remains after correcting for 
the parallactic motion and proper motion of our 
target. 


7.3. Astrometric Results 


In Table IH we list the number of seasons the 
target has been observed, the number of parallax 
frames, the start and end dates of observations, 
the time duration of the observations, the num¬ 
ber of reference stars, the relative parallax (tt 
rel), the correction to absolute parallax (tt corr), 
the absolute parallax (tt abs), the proper motion 
amplitude (/r), the proper motion position an¬ 
gle (P.A.), and the tangential velocity (vtan) for 
the 12 targets in our sub-sample. We do note 
that the parallax correction to absolute for GJ 
628 was much larger (4.8 (t) than typical. We ex¬ 
pect that this is a consequence of the reference 
stars being highly reddened, which will bias dis¬ 
tance estimates. Rather than using this number, 
we use the average of the corrections to abso- 


lute from 221 stars previously 

published using the 

CTIOPI nioeline (Jao et al. 

2005: Heur 

V et al. 

2006: Gizis et al. 2007: Subasavage et al. 

20091: 

Riedel et al. 2010l: l.lao et al.l 

201 ll: Riedel et al. 

2011 : Ivon Braun et al. 201l|). 


All of our measured parallaxes are within 3 
sigma of the weighted average of previously pub¬ 
lished values, except GJ 3( 30. Our paral l ax val ue of 
GJ 300 supersedes that in iHenrv et al.l ( 2006h . be¬ 
cause we now have roughly twice as many frames 
over twice the time span, and now use an improved 
centroiding technique. 


After solving for the parallactic and proper mo¬ 
tions, the average residual deviations for all 12 


stars is 2.74 mas in right ascension (R.A.) and is 
3.36 mas in declination (DecL). The average resid¬ 
ual deviations range from 0.86 mas to 4.85 mas in 
R.A. and from 1.54 mas to 5.53 mas in Decl. (see 
Table [T]) . Residual deviations are calculated by 
taking the standard deviation of the absolute val¬ 
ues of nightly mean positions. We also report the 
statistical uncertainty (henceforth referred to as 
the mean error) for each star, which ranges from 
2.39 mas to 6.18 mas in R.A. and from 2.65 mas 
to 6.88 mas in Decl. We calculate the error on a 
single night by taking the standard deviation of 
the offsets from zero for the frames, typically five, 
taken on each star during a night. Then, we cal¬ 
culate the mean error for the star by taking the 
average of all the nightly errors. The mean errors 
for both R.A. and Decl. are also listed Tabled On 
average our errors are 1.5 times larger than resid¬ 
ual deviations. We take this as indication that 
we are slightly over estimating our errors. For all 
12 stars, the astrometric signals that remain after 
correcting for the parallactic and proper motions 
are plotted over time in Figure 4, split into R.A. 
and Decl. 

8. Companion Detection Limits 

To identify companions around these 12 stars, 
we search for hidde: 
using Lomb-Scargle 
For each star, we search for companions with pe¬ 
riods between 2 to 100 days for our RV data and 
between 300 and 3000 days for our astrometric 
data (in R.A. and Decl. separately). Using the 
IDL program scargle.pro, no frequencies have pow¬ 
ers that exceed the Icr false alarm probability. We 
conclude that there are no periodic signals in ei¬ 
ther our RV or astrometric data. From this, we 
assume our stars are single stars within our detec¬ 
tion limits. 

To set limits on the presence of companions, we 
perform Monte Carlo simulations to determine the 
minimum object mass that we would have been 
able to detect in our data given our measurement 
errors and observing cadence. 

8.1. Radial Velocity Limits 

As reported in §6.4, these stars all have con¬ 
stant RVs to within 139 m s“^. To determine 
the minimum mass companion we would expect 


1 periodic sig nal in the dat a 
scarg 
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to detect, we simulate 1,000,000 circular orbits for 
each star allowing the inclinations and companion 
masses to vary at orbital periods between 0.5 and 
100 days. For each star, we use the stellar masses 
calculated in §2 and allow the stellar masses to 
fluctuate within 20% to account for possible er¬ 
rors in the mass estimates. We assume circular or¬ 
bits, which is likely appropriate for com panions to 
old st ars with short periods 10 davs: lLin fc Gu 
2 OO 4 II . The cosine of inclination ranges from 0 


to 90 degrees. After creating randomly orientated 
circular orbits, we then extract predicted radial ve¬ 
locity values at the Julian Dates we observed from 
these simulated orbits. Then, we compare these 
extracted radial velocity values to our determined 
RV errors at their respective Julian Date (see Ta¬ 
ble |4]) to determine a value. Then, we use 
this x^ valu e as an input for the IDL routine, mp- 
chitest.pro I Markwardtll2009l l. to determine which 
orbits we would have been able to detect with a 
3a confidence. Figures 5-16 provide illustrations 
of the RV Monte Carlo simulation results for our 
stars. The plots are shaded to display the fraction 
of objects that would have been detected with a 
3a confidence in each period and mass increment 
bins. The sizes of the mass and period bins are 
0.10 Mji/p and 1 day, respectively. The illus¬ 
trations show that most Jupiter-mass companions 
with periods less than 30 days would be detected. 
Although some Jupiter-mass companions could be 
detected at longer periods, the sensitivity becomes 
strongly dependent on the observing cadence rela¬ 
tive to the orbital ephemeris. Aliasing is apparent 
in our RV data. 


To interpret our data collectively, we follow the 
same method as above, except that we determine 
the companion mass that we would have been able 
to detect at selected orbital periods of 3 days, 10 
days, 30 days, and 100 days. We allow our se¬ 
lect period to range from ±10% to mitigate the 
bias caused by aliasing and give a better estimate 
of the masses of the systems we would be able to 
detect. The minimum companion mass we would 
have been able to detect is given in increments of 
0.5 Mjup. On average for our 12 stars, these sim¬ 
ulations suggest that with 3a confidence, we would 
have been able to detect 90% of companions with 
a mass of ~1 Mjpp in 3 day orbits, ~1.5 Mjpp 
in 10 day orbits, ~3 Mjpp in 30 day orbits and 
~7 Mjj/p in 100 day orbits. The individual detec¬ 


tion limits for each star are provided in Table [8] 
GJ 1065 is only observed over a baseline of ~13 
days, hence the limits on the mass of a companion 
that we would be able to detect around this star 
are significantly larger for an orbital period of 100 
days than the rest of stars in this sub-sample. 


Exoplanets with longer orbital periods can ex¬ 
hibit large eccentricities. In our orbital period 
range of 0.5 to 100 days, there are 1356 exoplanets 
known to dat Of those, close to 5% have ec¬ 
centricities listed above 0.2 and less than 1% of ex¬ 
oplanets have eccentricities above 0.6. The planet 
with highest eccent ricity in this period r ange is Ke¬ 
pler 419b (e=0.85; Dawson et al. 2014 ). To be in¬ 
clusive of planets with non-zero eccentricities that 
have orbital periods between 30 and 100 days, we 
performed our simulations again allowing the ec¬ 
centricity to vary from 0.0 to 0.9. In these cases, 
we find that the minimum companion mass we 
could detect with our data with a 90% detection 
rate increases by a factor of three when the eccen¬ 
tricity varies uniformly over this range, compared 
to when the eccentricities are set to zero. 


8.2. Astrometric Limits 


Similar to the Monte Carlo simulations per¬ 
formed for the RV orbits, we simulate 1,000,000 
photocentric orbits for each star allowing the in¬ 
clinations, companion masses and eccentricities to 
change. For each star, we use the stellar masses 
calculated in §2 and allow the stellar masses to 
fluctuate within 20%. We allow the inclination to 
vary between 0 to 90 degrees and the eccentricity 
to vary between 0 and 0.9. We allow the longitude 
of periastron, oj, and the longitude of the ascend¬ 
ing node, D, to vary from 0 to 180 degrees. We use 
the weighted distance measurements determined 
from the literature and from our parallax measure¬ 
ment in §7.3. We allowed the orbital periods to 
range from 2 to 8 years to represent our observing 
cadence and temporal baseline. We then use these 
simulated orbits to extract astrometric positions 
at our observed Julian Dates. We compare these 
extracted astrometric positions to our determined 
astrometric errors (see Table [7]) to determine a x^ 
value. Then, we use this x^ va lue as an input fo r 
the IDL routine, mpchitest.pro I Markwardt 200^, 
to determine which orbits we would have been able 


^'^Numbers are from exoplanet.eu as of October 27, 2014. 
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to detect with a Scr confidence. Figures 5-16 dis¬ 
play the astrometric Monte Carlo simulation re¬ 
sults. The plots are shaded to include the fraction 
of objects that would have been detected with a 
3 ct confidence in each period and mass increment 
bin of 0.10 Mjjjp and 0.1 years, respectively. 

To obtain a quantitative estimate of our de¬ 
tection limits, we again perform 1,000,000 Monte 
Carlo simulations. We keep the same parameters 
as described above, except we search for compan¬ 
ions at selected orbital periods of 2 years, 4 years, 
6 years, and 8 years. We give the minimum com¬ 
panion mass we would have been able to detect in 
increments of 0.5 Mjf/p. On average for our 12 
stars, these simulations suggest that with 3 (t con¬ 
fidence, we would have been able to detect 90% of 
companions with masses of ~18.5 Mjpp in 2 year 
orbits, ~11.5 Mjpp in 4 year orbits, ~9 Mjpp in 
6 year orbits and ~7.5 Mjpp in 8 year orbits for 
all of our stars (see Table [S]) . 


8.3. Additional Detection Limits from the 
Literature 

8.3.1. Optical RV Limits 


To place the strictest limits on excluding com¬ 
panions for the stars in this sub-sample, we also 
include the results for eleven of the stars that have 
been previously observed with high dispersion op¬ 
tical spectroscopy to sear ch for planetary compan¬ 
ions. iBarnes et al.l ( 2012ll present 8 epochs of GJ 
1286 with a root mean square scatter of 22.1 m 
s“^ over 1.1 d ays. With th a t prec ision and short 
time cadence, Barnes et al. ( 2012h would only be 
sensitive to planets with periods on the order of a 
few days. The remaining ten stars with previous 
RV measurem ents are included in the HARPS M- 
dwarf survey (|Bonfils et al.l 120131 1. Of those ten, 
the authors present Monte Carlo results to exclude 
the presence of companions around eight of the 
stars, as the other two stars have fewer than four 
epochs of data. On average for the eight stars, the 
authors are able to exclude with a 99% confidence 
the presence of companions with masses greater 
than 0.2 Mjpp in 10 to 100 day orbital periods 
and greater than 1.3 M jpp in an eight year or¬ 
bit (|Bonfils et al.ll2013^ . We find no previous high 
dispersion spectroscopy data for G 99-49; thus our 
data place the strictest companion limits for this 
star. 


8.3.2. High Resolution Imaging Limits 


To include companions beyond several AU, we 
include high resolution results for companions. 
AO imaging or HST imaging to search for stel¬ 
lar and brown companions has been completed 


2012 : 


2001 


for 11 out of the 12 stars (lDietericlL_et_aL 
NaJ«ajjrnn_et ah _ H194I: ODpenheimer et al.l 
Tanner et al ]l2om GJ 1065 has not been imaged 
with high resolution techniques . For 10 of the 12 
stars. lOppenheimer et al.l (|200lll would have been 
able to detect any stellar companions at separa¬ 
tions greater than 10 AU and most old (<5 Gyr) 
brown dwarfs with ma sses more than 40 M./g p be¬ 


tween 40 and 120 AU. Dieterich et al. (2012!) have 


further constrained the spatial and mass regime 
of comp anions around 7 of the 12 stars in the sub¬ 
sample. Dieterich et al.l ( 20121) adopt 0.04 Msun 
(~42 Mj( 7 p) at 3 Gyr as the minimum mass de¬ 
tectable with their HST/NICMOS snapshot high 
resolution imaging survey, and find no L dwarf 
companions in the separation range of 5 to 70 AU 
and no T dwarf companions (~0.05-0.08 Msj/at) 
to these M-dwarfs in the separation range of 10 to 
70 AU. 

More targeted surveys have pushed the de¬ 
tection limits of companions down to the plan- 
etary range with the use of near infrared AO 


j Lafreniere_e^^ ■l2007l:l iLuhman et al.ll2005l:lMasci adri etM 


20051: Schroeder et al. 1*20001 ). iMasciadri et ah 
(12005 ) find that GJ 628 has no companions more 
massive than 5 Mjup at 4 AU and no com- 
panions more massive than 10 Mjjjp at 1 AU. 
Lafreniere et al.l (j2007^ also rule out the presence 
of companions more massive than 5 Mjpp in the 
separation range of 25 to 50 AU around GJ 628. 
These two studies imply that we can also find mas¬ 
sive planets between 1 and 50 AU around nearby 
M-dwarfs via imaging techniques, assuming the 
luminosities from models are correct. 


8.3.3. Common Proper Motion Limits 

To extend our search to inlcude wide compan¬ 
ions (> 100 AU), we include co mmon proper mo 


tion searches for companions. iHinz et al.l (12002^ 


have searched 10 of the 12 targets and find no stel¬ 
lar or brown dwarf companions at wide separations 
(~100—1400 AU) around these stars. With a lim¬ 
i ting J band magnitude of ~16.5 mag. iHinz et al 
(2002) would have been able to find ~40 Mjjjp 
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for an assumed age of 5 Gyr at a distance of 5.8 
pc. 


9. Summary 


In recent years, we have gained a better statis¬ 
tical understanding on how common Jupiter-size 
and Earth-size planets are around stars from the 
Kepler mission that photometrically monitored 
appro ximately 150,000 sta rs in its first 3.5 years 
' e.g. Batalha et ^ 2010 : Howard et al. 2012t 


Dressing fc CharbonnearJ2013 : Fressin et al. 2013 


Mulders et al .11201411 . Of those, 3,000 are M-stars 
that are brighter than 1 6th magnitude in th e Ke¬ 


pler passband (Kp<16 Batalha et al. 2010ll . As 


Kepler is magnitude limited, most M-stars being 
monitored are early M-dwarfs. Kepler has dis- 
coyered many planets around these stars and even 
yound mid M-dwarfs, like t he M4 dwarf Kepler 42 
( Muirhead et al.l2012( 201^ . However Keyler and 


the subsequent K2 mission (jHowell et al.l 120141) 


will not able to survey all of our faint stellar neigh¬ 
bors and thereby determine an unbiased measure¬ 
ment of the frequency of planets around the latest 
M-dwarfs. With the current limited number of 
mid to late M-dwarfs that have been surveyed 
for companions by ground based telescopes and 
Kepler, we cannot concretely say whether planets 
even exist around these stars, as the latest known 
star with an exoplanet is the M4.5 dwarf G J 
1214 ( Reid et al. 1995 : Gharbonneau et al. 20091) . 


Therefore, the New Worlds New Horizons 2010 
Astronomy Decadal Survey has stated the need to 
develop new innovative reduction techniques and 
instrumentation to explore the lowest mass stars 
on the main sequence, as these stars are often 
too faint and chromospherically active for current 
ground based optical companion surveys and the 
Kepler mission. Motivated by this, we have iden¬ 
tified and characterized a volume-limited survey 
of mid to late M-dwarfs and present our infrared 
radial velocity and astrometric programs, which 
are both viable methods to discover Jupiter-mass 
and brown dwarf companions. 

At the start of this study, we identified a 
volume-limited equatorial sample of 60 mid M- 
dwarfs that extends out to 10 pc, which are observ¬ 
able from most facilities in the northern and south¬ 
ern hemispheres. Since this time, two of these 
stars have been found to have spectroscopic com- 


panions (GJ 867BD, LHS 1610AB: iDavison et al 


2014 : Bonfils et al.l 20131 ). For the remaining 58 


isolated stars in the sample, we provide new V, R 
and I photometry and new low dispersion opti¬ 
cal (6000 — 9000 A) spectra, from which we mea¬ 
sure the Ha equivalent widths, Nal indices and 
spectral types. We determine that 35 stars are 
emission-line stars, and note the trend of later 
type M-stars being more likely to be emission-line 
stars than earlier type M-stars. 


For a smaller sub-sample of 12 stars, we present 
rotational velocities and absolute radial velocities. 
Ten of these 12 stars have vsini values below our 
detection limit (3 km s“^), while G 99-49 and 
GJ 729 have vsini values of 5.8+0.3 km s“^ and 
3.8+0.6 km s“^, respectively. For these 12 stars 
our observed RV dispersion is 99 m s“^ with a 
standard deviation of 27 m s“^. We also demon¬ 
strate the achievable RV precision using our tech¬ 
nique is ~90 m s“^ over timescales from 13 days 
to 5 years. Our spectroscopic results indicate that 
on average we would have detected companions 
with masses greater than 1.5 Mjjjp in 10 day or¬ 
bital periods and greater than 7 Mjup in 100 day 
orbital periods. 


Over baselines of 9—13 years, we provide as¬ 
trometry with typical residuals (after determining 
parallaxes and proper motions) of ~3 mas. This 
allows us to exclude the presence of companions 
with masses greater than 11.5 Mjyp with orbital 
periods of 4 years and 7.5 Mjup with orbital pe¬ 
riods of 8 years. Although we do not detect com¬ 
panions around any of the stars in our sub-sample, 
these results do show that we could easily detect 
brown dwarfs with wide orbits around low mass 
stars. 


An ensemble result of this work is that all 12 
stars studied here with both infrared spectroscopy 
and astrometry and in the literature with various 
other techniques are found to be single stars with¬ 
out stellar, brown dwarf, or giant planet compan¬ 
ions within the respective detection limits of the 
studies. These results provide a first step in the 
process of vetting our nearest neighbors for planet 
searches. We suggest these slowly rotating, single 
stars would be prime targets for upcoming Earth 
mass RV planet searches using new instruments 
(e.g. ISHELL with use of a gas ce l l, SPIROU and 


CARMENES iR.a vner et al.ll2012t iReshetov et al 
2012 : Qnirrenbach et al .l l201(C coming online, as 
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these instruments should be sensitive enough to 
search for Earth-mass planets around mid to late 
M-dwarfs. 
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Table 1 

Effectively Single Equatorial (±30° Decl.) mid M-dwarfs within 10 pc 


Name 

R.A. J2000.0 

Decl. 

J2000.0 

Parallax^ 

(mas) 

SpType 

Mass^ 

(Mq) 

Mv 

(mag) 

Vj 

(mag) 

Rkc 

(mag) 

Ikc 

(mag) 

# nts. 

r 

(mag) 

H‘= 

(mag) 

(mag) 

(km 

GJ 1002 

00 

06 

43.2 

-07 

32 

17.0 

213.00+3.601 

M5.0V 

0.11 

15.48 

13.84 

12.21 

10.21 

3 

8.23 

7.79 

7.44 

<3.0i 

GJ 54.1 

01 

12 

30.6 

-16 

59 

56.3 

269.08+2.99i’2 

M4.0Ve 

0.13 

14.30 

12.15 

10.73 

8.95 

2 

7.26 

6.75 

6.42 

<2.5® 

LHS 1302 

01 

51 

04.5 

-06 

07 

04.8 

100.78+1.893 

M4.5Ve 

0.13 

14.51 

14.49 

13.00 

11.16 

5 

9.41 

8.84 

8.55 


GJ 83.1 

02 

00 

13.0 

+ 13 

03 

07.0 

224.80+2.901 

M4.0Ve 

0.14 

14.11 

12.35 

10.95 

9.18 

3 

7.51 

6.97 

6.65 

<2.5® 

LHS 1326 

02 

02 

16.2 

+ 10 

20 

13.7 

112.00+3.201 

M5.0Ve 

0.10 

15.95 

15.70 

13.99 

11.91 

3 

9.84 

9.25 

8.93 

<4.5® 

LHS 1375 

02 

16 

29.8 

+13 

35 

13.7 

117.70+4.001 

M5.0Ve 

0.10 

16.15 

15.80 

14.14 

12.01 

3 

9.87 

9.31 

8.98 

12.4® 

GJ 102 

02 

33 

37.2 

+24 

55 

39.2 

102.40+2.701 

M3.5Ve 

0.17 

13.10 

13.05 

11.76 

10.10 

3 

8.47 

7.91 

7.63 


GJ 105B 

02 

36 

15.4 

+06 

52 

19.1 

138.79±0.43i’2A 

M3.5V 

0.22 

12.42 

11.71 

10.49 

8.88 

2 

7.33 

6.79 

6.57 

<2.5® 

SO 0253+1652 

02 

53 

00.8 

+ 16 

52 

53.3 

259.41+0.893’5 

M7.0V 

0.08 

17.21 

15.14 

13.03 

10.65 

3 

8.39 

7.88 

7.59 

10.0® 

LP 771-095A 

03 

01 

51.1 

-16 

35 

30.7 

143.81+2.493 

M2.5V 

0.25 

12.01 

11.22 

10.07 

8.66 

4 

7.11 

6.56 

6.29 

5.5® 

GJ 1057 

03 

13 

23.0 

+04 

46 

29.3 

117.10+3.501 

M4.5V 

0.13 

14.28 

13.94 

12.45 

10.62 

3 

8.78 

8.21 

7.83 

<2.2® 

GJ 1065 

03 

50 

44.3 

-06 

05 

41.7 

105.40+3.201 

M3.5V 

0.18 

12.93 

12.82 

11.60 

10.04 

3 

8.57 

8.00 

7.75 

4.0® 

GJ 166G 

04 

15 

21.7 

-07 

39 

17.4 

200.65+0.23^ 

M4.0Ve 

0.19 

12.75 

11.24 

9.99 

8.31 

3 

6.75 

6.28 

5.96 

5.0® 

LP 655-048 

04 

40 

23.3 

-05 

30 

08.3 

105.50+3.208 

M6.0Ve 

0.08 

17.92 

17.80 

15.73 

13.37 

5 

10.66 

9.99 

9.55 

16.5® 

LHS 1723 

05 

01 

57.4 

-06 

56 

46.5 

187.92+1.263 

M4.0Ve 

0.15 

13.57 

12.20 

10.86 

9.18 

4 

7.62 

7.07 

6.74 

<3.2®" 

GJ 203 

05 

28 

00.2 

+09 

38 

38.1 

102.60±2.09i’2 

M3.0V 

0.19 

12.52 

12.46 

11.27 

9.78 

3 

8.31 

7.84 

7.54 

4.0® 

GJ 213 

05 

42 

09.3 

+ 12 

29 

21.6 

171.50+1.OOi’^T 

M3.5V 

0.19 

12.71 

11.54 

10.32 

8.68 

2 

7.12 

6.63 

6.39 

<2.5® 

G 99-49 

06 

00 

03.5 

+02 

42 

23.6 

190.77+1.86i’3 

M3.5Ve 

0.19 

12.71 

11.31 

10.04 

8.43 

6 

6.91 

6.31 

6.04 

7.4'® 

GJ 232 

06 

24 

41.3 

+23 

25 

58.6 

119.40+2.301 

M4.0V 

0.15 

13.55 

13.16 

11.86 

10.21 

3 

8.66 

8.16 

7.91 

<3.1® 

GJ 1093 

06 

59 

28.7 

+ 19 

20 

57.7 

128.80+3.501 

M5.0Ve 

0.11 

15.49 

14.94 

13.25 

11.24 

4 

9.16 

8.55 

8.23 

<2.8® 

GJ 273 

07 

27 

24.5 

+05 

13 

32.8 

266.23+0.66i’2.7 

M3.0V 

0.25 

12.01 

9.88 

8.68 

7.14 

3 

5.71 

5.22 

4.86 

2.5® 

GJ 283B 

07 

40 

19.2 

-17 

24 

45.0 

109.45+0.511’®’® 

M6.5Ve 

0.16 

13.26 

13.06 

12.89 

12.72 

4 

10.16 

9.63 

9.29 


GJ 285 

07 

44 

40.2 

+03 

33 

08.8 

167.19+2.05i’2 

M4.0Ve 

0.23 

12.31 

11.19 

9.91 

8.22 

4 

6.58 

6.01 

5.70 

4.5®® 

GJ 1103 

07 

51 

54.7 

-00 

00 

11.8 

114.00+3.301 

M4.5V 

0.15 

13.54 

13.26 

11.89 

10.19 

3 

8.50 

7.94 

7.66 


GJ 299 

08 

11 

57.6 

+08 

46 

22.1 

146.30+3.101 

M3.5V 

0.15 

13.69 

12.86 

11.57 

9.91 

3 

8.42 

7.93 

7.66 

3.0'^ 

GJ 300 

08 

12 

40.9 

-21 

33 

06.8 

125.78+0.973 

M3.5V 

0.20 

12.65 

12.15 

10.85 

9.22 

3 

7.60 

6.96 

6.71 

<3.0® 

GJ nil 

08 

29 

49.3 

+26 

46 

33.7 

275.80+3.001 

M6.0Ve 

0.09 

17.16 

14.96 

12.89 

10.59 

3 

8.24 

7.62 

7.26 

8.E 

LHS 2090 

09 

00 

23.6 

+21 

50 

05.4 

156.87+2.673 

M6.0Ve 

0.09 

17.09 

16.11 

14.12 

11.84 

3 

9.44 

8.84 

8.44 

20.0® 

LHS 2206 

09 

53 

55.2 

+20 

56 

46.0 

108.69+2.063’10 

M4.0Ve 

0.14 

14.20 

14.02 

12.63 

10.85 

3 

9.21 

8.60 

8.33 

16.5® 

LHS 292 

10 

48 

12.6 

-11 

20 

08.2 

220.30+3.601 

M6.5Ve 

0.08 

17.50 

15.78 

13.63 

11.25 

3 

8.86 

8.26 

7.93 

3.0® 

GJ 402 

10 

50 

52.0 

+06 

48 

29.2 

145.67+3.17i’2 

M4.0V 

0.21 

12.53 

11.71 

10.43 

8.84 

3 

7.32 

6.71 

6.37 

<2.5® 

GJ 406 

10 

56 

28.9 

+07 

00 

53.2 

419.10+2.101 

M5.0Ve 

0.09 

16.69 

13.58 

11.64 

9.44 

2 

7.09 

6.48 

6.08 

<3.0® 

GJ 447 

11 

47 

44.4 

+00 

48 

16.4 

298.14+1.37i’2 

M4.0V 

0.16 

13.52 

11.15 

9.79 

8.13 

3 

6.51 

5.95 

5.65 

<2.5® 

GJ 1154 

12 

14 

16.5 

+00 

37 

26.4 

119.40+3.501 

M4.5Ve 

0.14 

14.03 

13.64 

12.17 

10.31 

2 

8.46 

7.86 

7.54 

5.2® 

GJ 1156 

12 

18 

59.4 

+ 11 

07 

33.9 

152.90+3.001 

M4.5Ve 

0.12 

14.87 

13.95 

12.33 

10.38 

2 

8.53 

7.88 

7.57 

9.2'® 

GJ 486 

12 

47 

56.6 

+09 

45 

05.0 

119.58+2.64® 

M4.0V 

0.28 

11.79 

11.40 

10.21 

8.67 

2 

7.20 

6.67 

6.36 

<2.5® 

GJ 493.1 

13 

00 

33.5 

+05 

41 

08.1 

123.10+3.501 

M4.5Ve 

0.15 

13.85 

13.40 

12.02 

10.26 

2 

8.55 

7.97 

7.66 

16.8'® 

GJ 555 

14 

34 

16.8 

-12 

31 

10.3 

160.78+1.98i’®’ii 

M3.5V 

0.22 

12.37 

11.34 

10.06 

8.44 

3 

6.84 

6.26 

5.94 

<2.5® 






Table 1 —Continued 


Name 

R.A. J2000.0 

Decl. 

J2000.0 

Parallax" 

(mas) 

SpType 

Mass*’ 

(M©) 

Mv- 

(mag) 

Vj 

(mag) 

Rkc 

(mag) 

Ikc 

(mag) 

nts. 

(mag) 

(mag) 

(mag) 

usim^ 
(km s“^) 

LHS 3003 

14 

56 

38.3 

-28 

09 

47.4 

152.494:2.02l’®'l2 

M7.0Ve 

0.08 

17.99 

17.07 

14.92 

12.54 

5 

9.97 

9.32 

8.93 

5 . 0 ® 

GJ 609 

16 

02 

51.0 

-120 

35 

21.8 

100.304:3.101 

M4.0V 

0.20 

12.59 

12.58 

11.32 

9.70 

2 

8.13 

7.65 

7.37 

< 3 . 0 ’*’ 

GJ 628 

16 

30 

18.1 

-12 

39 

45.4 

234.384:1.50i’2 

M3.0V 

0.26 

11.92 

10.07 

8.89 

7.37 

3 

5.95 

5.37 

5.08 

1.510 

GJ 643 

16 

55 

25.3 

-08 

19 

20.8 

154.964:0. 52 L 8 . 13 ,14 

M3.0V 

0.19 

12.72 

11.77 

10.55 

9.01 

3 

7.56 

7.06 

6.72 

< 2 . 7 ’*' 

GJ 644G 

16 

55 

35.3 

-08 

23 

40.1 

154.964:0. 52 L 8 . 13 ,14 

M7.0Ve 

0.08 

17.80 

16.85 

14.64 

12.25 

3 

9.78 

9.20 

8.82 

9.01 

GJ 1207 

16 

57 

05.7 

-04 

20 

56.0 

115.264:1.503 

M3.5Ve 

0.19 

12.56 

12.25 

10.99 

9.43 

5 

7.97 

7.44 

7.12 

10 . 7 ® 

GJ 699 

17 

57 

48.5 

4-04 

41 

36.2 

545.514:0.29i’^'i® 

M3.5V 

0.17 

13.17 

9.49 

8.27 

6.70 

1 

5.24 

4.83 

4.52 

<2.5^ 

GJ 1224 

18 

07 

32.9 

-15 

57 

47.0 

132.604:3.701 

M4.0Ve 

0.14 

14.09 

13.48 

12.08 

10.31 

3 

8.64 

8.09 

7.83 

<3.0io 

GJ 1230B 

18 

41 

09.8 

4-24 

47 

19.5 

120.904:7.201 

M4.5Ve 

0.19 

12.74 

12.33 

10.97 

9.26 

4 

8.86 

8.0 

7.77 

<7.li 

GJ 729 

18 

49 

49.4 

-23 

50 

10.4 

337.224:1.97i’2 

M3.5Ve 

0.17 

13.14 

10.50 

9.26 

7.68 

3 

6.22 

5.66 

5.37 

4.0= 

GJ 752B 

19 

16 

57.6 

4-05 

09 

02.2 

171.204:0.SQi’^.ie,17 

M8.0Ve 

0.07 

18.62 

17.45 

15.21 

12.78 

4 

9.91 

9.23 

8.77 

6 . 5 I 

GJ 1235 

19 

21 

38.7 

4-20 

52 

02.8 

100.104:3.501 

M4.0Ve 

0.16 

13.47 

13.47 

12.12 

10.46 

2 

8.80 

8.22 

7.94 


GJ 1256 

20 

40 

33.6 

-115 

29 

57.2 

102.004:2.201 

M4.0V 

0.16 

13.51 

13.47 

12.10 

10.37 

3 

8.64 

8.08 

7.75 

< 6.511 

LP 816-060 

20 

52 

33.0 

-16 

58 

29.0 

175.034:3.402 

M3.0V 

0.19 

12.72 

11.50 

10.25 

8.64 

3 

7.09 

6.52 

6.20 

< 6.511 

G 188-038 

22 

01 

13.1 

4-28 

18 

24.9 

111.704:1.73i’2 

M3.5Ve 

0.23 

12.29 

12.05 

10.77 

9.16 

2 

7.64 

7.04 

6.78 

35.1® 

LHS 3799 

22 

23 

07.0 

-17 

36 

26.1 

134.404:4.901 

M4.5Ve 

0.14 

13.94 

13.30 

11.87 

10.04 

5 

8.24 

7.64 

7.32 

< 6.511 

LP 876-010 

22 

48 

04.5 

-24 

22 

07.5 

132.074:1.191® 

M4.0Ve 

0.17 

13.19 

12.59 

11.31 

9.61 

3 

8.08 

7.53 

7.21 

22.012 

GJ 896A 

23 

31 

52.2 

-119 

56 

14.3 

159.884:1.53i’2.i9 

M3.5Ve 

0.33 

11.32 

10.30 

9.13 

7.66 

2 

6.16 

5.57 

5.33 

10.0® 

GJ 896B 

23 

31 

52.6 

-119 

56 

13.9 

159.884:1.53i’2.i9 

M4.0Ve 

0.16 

13.42 

12.40 

11.04 

9.28 

2 

7.10 

6.56 

6.26 

15.0® 

GJ 1286 

23 

35 

10.5 

-02 

23 

20.8 

138.304:3.501 

MS.OVe 

0.11 

15.43 

14.73 

13.10 

11.10 

3 

9.15 

8.51 

8.18 

< 5.71 


^When mu l tiple p aral lax references are li sted, the reported value here is the weighted means for e ach s ystem. Parallax Refer ence s: Ivan A ltena et al.l lll995[l. (2) 


When mu l tiple p aral lax reier 

Ivan LeeuwenI |[2007tl: ('3)lHenrv et al.l ||200€ h f4)llan na" et aL| ljl996l ') ; CS) iGate wood fc Cobar _j2009|h (6) fShkolnik et ^ 

1I2OO5I): ^9) ISub asavag e et al. | l|2009l'). flO ) Smart et alj ||201C|): fll) iJao et al.| 1I2OO5I) ; (liy Tinned lil996h: (13) ISoderhieln ^ _ , ^ ^ ___ „ _ ; , , 

iBenedict et al.l l|l999l V: G6) Ixinnev et al] 1 1995h : ( 17^ IPravdo &: Shaklanl ll20odh (18) lMamaiek et al.l l|2013fl : G9) IWeid lIlQQa '). This table only includes previously 
published parallax measurements. New astrometric measurements from this paper are given in Table 6. 


h : ('7') | Gatewoodl ||2008|): . 

199^: (14‘1 iMartin et al.l 


ICos^ et al.l 

199^); (15) 


^M asses were determined using the mass-luminosity relationship of lHenrv et al.l l|l999tl for stars with > 12.89 mag, and the relationship of lHenrv McCarthvI 
{TUI) for brighter stars. Typically errors using these relations are close to 20%. 

^All J, Ks magnitudes are from the 2MASS All Sky Catalogue of point sources from iskrutskie et al.l <l2006h . 


^usim R eferences: fl^lMohantv &; Bas Til 1120031): ('2') iBrowning et al.l <l20 10h: f 3 ^|Jenkins et al.l <l2009lh {'4)lBarnes et al.l ll20l4): fS^lTanner et all ll2012lh (6) iReinersI 
2^R); (7) iDelfosse et ahl jllOOSlb f8) iJones et alj il2005l' ): f9) iReiners &; BasrJ ll2010ri : TO') iReiners &; Basril ll2007f) : Til) iBonfils et al.l ll2013lb (12) iMamaiek et al. 


2013). 




















































































Table 2 

Close-separation (<T0 M-Dwarf Equatorial (±30° Decl.) Multiples within 10 pc 


Name 

R.A. J2000.0 

Decl. J2000.0 

Parallax^ 

(mas) 

SpType® 

(mag) 

yd 

(mag) 

R'^ 

(mag) 

jd 

(mag) 

J'* 

(mag) 

Rd 

(mag) 

As'! 

(mag) 

Configuration^ 

GJ 1005AB 

00 15 27.7 

-16 07 56.0 

168.42+0.89i’2.3.4 

M4.0VJ1 

13.75 

12.621 

11.461 

IO.O5I 

7.22 

6.71 

6.39 

AB 

GJ 2005ABC 

00 24 44.2 

-27 08 25.2 

129.71±2.43i’® 

M5.5VJ2 

15.98 

15.42^ 

13.712 

11.562 

9.25 

8.55 

8.24 

ABC 

GJ 65AB 

01 39 01.5 

-17 57 01.8 

373.70±2.70i 

M5.5V1 

14.92 

12.06^ 

10.402 

8.342 

6.28 

5.69 

5.34 

AB 

GJ 105C 

02 36 04.7 

+06 53 14.8 

138.79+0.431’®-’^ 

M7.0V® 








AC-B 

LP yTi-ogeBC-f 

03 01 51.4 

-16 35 36.1 

143.81+2.49®’® 

M3.50VJ'‘ 

12.16 

11.37® 

10.13® 

8.58® 

7.29 

6.77 

6.50 

A-BC^ 

LHS 1610AB 

03 52 41.7 

+17 01 05.7 

100.88+2.05® 

M4.0VJ‘i 

13.87 

13.85® 

12.42® 

10.66® 

8.93 

8.38 

8.05 

AB 

GJ 190AB 

05 08 35.1 

-18 10 19.4 

107.57+2.081’® 

M3.5VJ2 

10.48 

10.32'! 

9 . 17 ! 

7 . 67 ! 

6.17 

5.59 

5.31 

AB 

GJ 234AB 

06 29 23.4 

-02 48 50.3 

244.44+0.92l’2’9 

M4.5VJ‘i 

13.06 

11.12'! 

9.78! 

8.08! 

6.38 

5.75 

5.49 

AB 

LTT 17993AB 

07 36 25.1 

+07 04 43.2 

116.60+0.97® 

M4.5VJ2 

13.58 

13.25® 

11.81® 

9.97® 

8.18 

7.61 

7.28 

AB 

GJ 1116AB 

08 58 15.2 

+19 45 47.1 

191.20+2.501 

M5.5VJ1 

15.06 

13.65® 

11.97® 

9.83® 

7.79 

7.24 

6.89 

AB 

LTT 12352AC 

08 58 56.3 

+08 28 25.9 

147.66+1.98® 

M3.5VJ‘i 

11.77 

10.92® 

9.67® 

8.05® 

6.51 

5.97 

5.69 

AC-B 

GJ 473AB 

12 33 16.3 

+09 01 26.0 

227.90+4.601 

M5.5VJ1 

14.28 

12.49'! 

10 . 93 ! 

8 . 97 ! 

7.00 

6.40 

6.04 

AB 

GJ 569BC 

14 54 29.4 

+16 06 08.9 

100.62+1.281’® 

M8.5VJ5 








A-BC 

GJ 644ABD 

16 55 29.6 

-08 19 55.3 

154.96+0.52i’2.5’io 

M2.5VJ1 

9.98 

9 . 03 ! 

7 . 94 ! 

6 . 57 ! 

5.27 

4.78 

4.40 

ABD-C-GJ 643 

GJ 695BC 

17 46 25.1 

+27 43 01.4 

120.32+0.161’® 

M3.5VJ2 

10.26 

9.86® 

8.70® 

7.25® 

5.77 

5.17 

4.95 

AD-BC 

GJ 1230AC 

18 41 09.8 

+24 47 14.4 

120.90+7.201 

M4.5VJ1 

12.57 

12.16® 

10.82® 

9.07® 

7.53 

6.91 

6.62 

AC-B 

GJ 791.2AB 

20 29 48.3 

+09 41 20.2 

112.90+0.301’® 

M4.5VJ® 

13.34 

13.08! 

11 . 73 ! 

9.98! 

8.23 

7.67 

7.31 

AB 

GJ 829AB 

21 29 36.8 

+17 38 35.9 

I 49 .OI+I. 69 I’® 

M3.4VJ1 

11.17 

10.30’^ 

9.15'^ 

7.70’^ 

6.25 

5.74 

5.45 

AB 

GJ 831AB 

21 31 18.6 

-09 47 26.5 

128.21+2. 05 I’® 

M4.5VJ1 

12.58 

12 . 04 ! 

10 . 74 ! 

9 . 04 ! 

7.32 

6.70 

6.38 

AB 

GJ 866ABC 

22 38 33.7 

-15 17 57.3 

289.50+4.401 

M5.0VJ1 

14.68 

12.372 

10.702 

8.642 

6.55 

5.95 

5.54 

ABC 

GJ 867BD 

22 38 45.3 

-20 36 51.9 

113 . 37 + 1 . 04 !’®’!! 

M3.5VJ1 

11.75 

11 . 45 ! 

10 . 29 ! 

8.78! 

7.34 

6.82 

6.49 

AC-BD 


When multiple parallax references are listed, the reported value here is the weighted mean for each system. Parallax References: (1) R^n ALterm_et_ajJ 


Soderhjehn|j [l^£ ); (3) ISmart et al.l _ _ _ _ , __ , 

Gatewood et alF 20031) : ('lO'l iMartin et al.1 l ll99al : ('iT fPavi^n et al. ||20U ). 


( 2 ) 


('4) [ Hershey Sz Tafj ||l99^ ) ; (5) Costa et aU ll2005l') : f6) lvan LeeuwenI ll2007j) : (7) llanna et al.l lll99dh fS) iHenT^^^jul ll200d) : (9) 

f ectral t 

ioii). 


*^9 represents a join t spectral type when two or m ore stars cannot be deco nvolved into their i ndivid ual components. Spectral type References: (1) iHenrv et aLl lll994|j : (2) 
iReid et al.l | |199'3) : ('3') iGolimowski et al.l ||2000|) : ('4‘): lHenrv et al] ll2006l) : fSj iKirkpatrick et al.l lll99ll) : fhj iDavison et ^ 


‘^Using our photometric measurements and the parallax data, the absolute magnitude errors range from 0.03 to 0.08 mag. 

“^We give joint photometry of the close binary stars as these close s ystems cannot be d econ volved into th eir in dividual componen ts. All stars with spect ral type den oted 
by the letter J also h ave jo int photometry. U, R and I References: flV lRiedel et al.l ll2010lj : f2j iBesse Iil lll99llj : f3j IHenrv et~^ ll200db ('4) [Besseil lll990lj : f5j IWeid lll996lj : (6) 
this paper; ('T) IWeij lll99lf) . All J, H, Ks magnitudes are from the 2MASS All Sky Catalogue of point sources from ISkrutskie et al.l ||2006 ). 


^System Configuration. We indicate widely (>4") separated pairs with a hyphen. There is no spacing between components with separations less than 4" (e.g. AB). 
^LP 771-096BC is a distant companion to LP 771-095, which is labeled A under this configuration. 

























































































Table 3 


Optical Spectroscopic Measurements for Effectively Single Equatorial mid M-dwarfs 


WITHIN 10 PC 


Name 

SpType 

Date 

YYYYMMDD 

H-a 
EW A 

Nal 

Index 

GJ 1002 

M5.0V 

19931031 

-0.1 

1.31 



20041001 

0.2 

1.36 

GJ 54.1 

M4.0Ve 

19931031 

-1.3 

1.25 



20041002 

-1.1 

1.30 

LHS 1302 

M4.5Ve 

20031206 

-3.6 

1.27 



20061209 

-3.4 

1.25 

GJ 83.1 

M4.0Ve 

19911113 

0.0 

1.20 



20041001 

-2.1 

1.28 

LHS 1326 

M5.0Ve 

20031206 

-1.4 

1.33 



20061208 

-0.8 

1.32 

LHS 1375 

M5.0Ve 

20031205 

-4.1 

1.33 



20061208 

-2.9 

1.32 

GJ 102 

M3.5Ve 

20031206 

-3.1 

1.22 

GJ 105B 

M3.5V 

19931029 

0.3 

1.12 



20041002 

0.2 

1.18 

SO 0253+1652 

M7.0V 

20031205 

0.8 

1.37 



20061206 

-0.4 

1.39 

LP 771-095A 

M2.5V 

20031206 

0.0 

1.16 

GJ 1057 

M4.5V 

20031206 

-0.3 

1.26 

GJ 1065 

M3.5V 

20020402 

0.0 

1.09 



20061207 

0.0 

1.21 

GJ 166C 

M4.0Ve 

20031206 

-5.3 

1.21 

LP 655-048 

M6.0Ve 

20040311 

-17.8 

1.35 

LHS 1723 

M4.0Ve 

20031206 

-1.5 

1.25 

GJ 203 

M3.0V 

19980208 

0.2 

1.13 



20061208 

0.1 

1.15 

GJ 213 

M3.5V 

19911113 

0.3 

1.12 



20040312 

0.0 

1.14 

G 99-49 

M3.5Ve 

19930314 

-3.4 

1.20 



20031206 

-4.6 

1.22 

GJ 232 

M4.0V 

19900122 

0.5 

1.17 



20040312 

0.0 

1.22 

GJ 1093 

M5.0Ve 

19930314 

-0.9 

1.29 



20031206 

-2.1 

1.33 

GJ 273 

M3.0V 

19900122 

0.4 

1.11 



20040312 

0.1 

1.15 

GJ 283B 

M6.5Ve 

19930316 

1.7 

1.33 



20031205 

-0.5 

1.40 



20061209 

-0.1 

1.38 

GJ 285 

M4.0Ve 

19930315 

-9.4 

1.20 



20041003 

-7.7 

1.24 

GJ 1103 

M4.5V 

20040312 

0.2 

1.20 

GJ 299 

M3.5V 

19930315 

1.6 

1.19 



20050131 

0.1 

1.24 

GJ 300 

M3.5V 

19930315 

0.8 

1.20 



20041003 

-0.2 

1.25 

GJ nil 

M6.0Ve 

19951202 

-4.4 

1.40 



20031206 

-3.8 

1.37 



20090203 

-4.0 

1.50 



20100301 

-7.6 

1.43 
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Table 3 —Continued 


Name 

SpType 

Date 

H-o 

Nal 



YYYYMMDD 

EW A 

Index 

GJ nil 

M6.0Ve 

20100305 

-8.0 

1.47 

LHS 2090 

M6.0Ve 

20020402 

-15.8 

1.16 



20061209 

-7.9 

1.33 

LHS 2206 

M4.0Ve 

20020401 

-4.4 

1.14 



20061209 

-3.3 

1.28 

LHS 292 

M6.5Ve 

19901123 

-0.8 

1.37 



20031207 

-2.4 

1.43 



20060531 

-9.5 

1.36 

GJ 402 

M3.0V 

20090203 

-0.2 

1.26 



20100305 


1.28 



20110517 

0.3 

1.21 

GJ 406 

M5.0Ve 

19951203 

-16.1 

1.37 



20030717 

-11.6 

1.33 



20090505 

-8.0 

1.38 

GJ 447 

M4.0V 

19930315 

1.0 

1.21 



20040608 

0.2 

1.23 



20090505 

-0.1 

1.26 

GJ 1154 

M4.5Ve 

20040313 

-6.6 

1.29 

GJ 1156 

M4.5Ve 

19930314 

-5.6 

1.25 



20040608 

-6.5 

1.27 

GJ 486 

M4.0V 

20040312 

0.3 

1.16 

GJ 493.1 

M4.5Ve 

20040312 

-4.9 

1.27 

GJ 555 

M3.5V 

19930314 

0.6 

1.18 



20040608 

0.5 

1.20 

LHS 3003 

M7.0Ve 

19930317 

-28.5 

1.29 



20040311 

-7.5 

1.40 



20060527 

-19.4 

1.35 



20090505 

0.0 

1.32 



20110727 

2.4 

1.00 

GJ 609 

M4.0V 

20090506 

-0.2 

1.16 



20110727 

0.1 

1.16 

GJ 628 

M3.0V 

19930315 

0.9 

1.13 



20040608 

0.0 

1.17 

GJ 643 

M3.0V 

20040608 

0.0 

1.19 

GJ 644C 

M7.0Ve 

19950812 

-11.6 

1.44 



20060525 

-7.1 

1.42 



20060526 

-10.4 

1.40 

GJ 1207 

M3.5Ve 

20020402 

-3.9 

1.11 



20060526 

-4.9 

1.17 



20090505 

-7.2 

1.27 

GJ 699 

M3.5V 

19950814 

0.0 

1.22 



20041001 

0.0 

1.21 



20060527 

0.6 

1.17 

GJ 1224 

M4.0Ve 

19931101 

-4.6 

1.25 



20040608 

-4.3 

1.26 

GJ 1230B 

M4.5Ve 

19931101 

-0.1 

1.26 



20040930 

-1.2 

1.29 

GJ 729 

M3.5Ve 

19931101 

-2.4 

1.20 

GJ 729 

M3.5Ve 

20040929 

-2.5 

1.25 

GJ 752B 

MS.OVe 

19950812 

-2.6 

1.31 
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Table 3 —Continued 


Name 

SpType 

Date 

YYYYMMDD 

H-a 
EW A 

Nal 

Index 



20030715 

-6.1 

1.33 



20040930 

-6.2 

1.29 



20060525 

-6.9 

1.27 

GJ 1235 

M4.0Ve 

20031011 

-0.5 

1.24 



20060526 

0.4 

1.21 

GJ 1256 

M4.0V 

20031011 

-0.2 

1.27 

LP 816-060 

M3.0V 

20020401 

0.5 

1.05 



20031207 

0.1 

1.19 

G 188-038 

M3.5Ve 

20031012 

-5.8 


LHS 3799 

M4.5Ve 

20031207 

-3.7 

1.30 



20060525 

-3.8 

1.31 



20060531 

-4.3 

1.27 

LP 876-010 

M4.0Ve 

20041002 

-6.5 

1.24 



20060525 


1.21 



20090726 

-3.3 

1.23 

GJ 896A 

M3.5Ve 

19931031 

-6.2 

1.17 



20041002 

-5.4 

1.17 



20110517 

-6.9 

1.17 

GJ 896B 

M4.0Ve 

19931031 





20110517 

-6.6 

1.17 

GJ 1286 

M5.0Ve 

19931031 

-0.7 

1.32 



20041001 

-0.7 

1.35 
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Table 4 

Radial Velocity Measurements 


Name 

HJD - 2,400,000 

Radial Velocity 

m 

SNR 

G 99-049 

54787.56* 

30171+ 92 

200 


54788.44 

30077+ 96 

258 


54790.51* 

30171+ 87 

254 


54791.47* 

30132+ 118 

202 


55146.52 

30244+ 122 

119 


55149.50 

30352+ 105 

140 


55151.53 

30317+ 102 

215 


55154.56 

30057+ 90 

208 


55638.31 

30144+ 95 

176 


55641.23 

30130+ 104 

145 


55642.27 

30294+ 86 

245 

GJ 300 

54790.60* 

9064+ 84 

199 


54791.65* 

9198+ 97 

103 


54963.27 

8882+ 85 

178 


54964.29 

9237+ 109 

89 


55151.64 

8922+ 96 

118 


55154.63 

9061+ 83 

184 


55317.25 

9047+ 83 

190 


55321.30 

8994+ 81 

181 


55644.36 

9042+ 84 

176 


56671.60 

8981+ 85 

154 

GJ 406 

54962.33 

19429+ 81 

212 


54963.38 

19415+ 83 

193 


54964.38 

19332+ 89 

148 


55320.41 

19759+ 81 

199 


55321.40 

19653+ 80 

204 


55323.41 

19548+ 84 

166 


55330.33 

19462+ 94 

87 


55331.30 

19465+ 80 

239 


55335.34 

19546+ 80 

162 


55637.46 

19731+ 80 

257 


55643.32 

19625+ 79 

180 


55645.39 

19489+ 80 

273 

GJ 555 

54959.43 

-1417+ 84 

214 


54960.44 

-1364+ 98 

116 


54961.41 

-1365+ 82 

197 


54962.44 

-1433+ 98 

238 


54963.47 

-1481+ 99 

124 


54964.50 

-1376+ 84 

166 


55320.52 

-1467+ 81 

199 


55355.44 

-1389+ 78 

226 

GJ 628 

54959.40 

-20969+ 89 

139 


54960.46 

-21135+ 96 

132 


54962.46 

-21058+ 80 

168 


54963.47 

-21100+ 89 

181 


54964.51 

-20980+ 91 

149 


55320.53 

-21257+ 87 

178 


55331.36 

-21205+ 77 

258 


55335.45 

-21180+ 83 

180 


56174.24 

-21095+ 79 

229 
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Table 4 —Continued 


Name 

HJD - 2,400,000 

Radial Velocity 
m s”^ 

SNR 

GJ 729 

54959.56 

-10109+ 84 

242 


54960.55 

-10178+ 82 

257 


54961.54 

-10236+ 97 

198 


54962.63 

-10127+ 109 

275 


55335.51 

-10419+ 81 

274 


55392.46* 

-10397+ 81 

228 

GJ 1002 

54787.41* 

-39770+ 85 

132 


54788.35 

-39886+ 101 

162 


54790.37* 

-39768+ 90 



55146.44 

-39912+ 95 

177 


55150.38 

-39976+ 86 

113 


55151.39 

-39913+ 84 


GJ 1065 

56259.53* 

-9118+ 87 

162 


56260.55 

-8984+ 92 

135 


56270.46 

-9163+ 95 

145 


56271.43 

-9122+ 137 

82 


56272.45 

-9110+ 93 

151 

GJ 1224 

54961.51 

-32425+ 127 

84 


54962.52 

-32609+ 135 

122 


54963.52 

-32811+ 96 

163 


54964.59 

-32559+ 164 

100 


55320.56 

-32560+ 90 

157 


55321.52 

-32704+ 90 

160 


55335.49 

-32789+ 87 

175 

GJ 1286 

54787.37* 

-40687+ 125 

74 


55146.37 

-40871+ 100 

126 


55149.29 

-40897+ 126 

76 


55151.36 

-40850+ 98 

109 


55154.40 

-40714+ 100 

138 

LHS 1723 

54787.51* 

42533+ 94 



54788.40 

42559+ 93 

149 


54790.47* 

42497+ 987 

188 


55149.47 

42327+ 90 

157 


55151.50 

42377+ 86 



55154.49 

42287+ 82 

217 

LHS 3799 

54790.28* 

-1643+ 104 

158 


55146.32 

-1568+ 92 

129 


55149.24 

-1612+ 83 

188 


55150.27 

-1651+ 85 

173 


55151.27 

-1596+ 90 

143 


55392.58 

-1779+ 88 

193 


*No telluric standards were observed on this night. 
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Table 5 

Spectroscopic Results 


Name 

Model Teff 
(K) 

Abs. RV“ 
(m s~^) 

N 

(# nights) 

A time 
(days) 

(7RV 
(m s“^) 

usim^ 

(km s”^) 

G 99-49 

3200 

30190 

11 

855 

98 

5.8+0.3 

GJ 300 

3200 

9043 

10 

1881 

110 

O.9+O.7I 

GJ 406 

3000 

19537 

12 

638 

132 

I.6-I-O.7I 

GJ 555 

3200 

-1413 

8 

396 

47 

O.6+O.7I 

GJ 628 

3000 

-21109 

9 

1215 

98 

I.2+O.7I 

GJ 729 

3200 

-10244 

6 

433 

125 

3.8+0.6 

GJ 1002 

3000 

-39871 

6 

368 

84 

I.4+O.3I 

GJ 1065 

3200 

-9099 

5 

13 

68 

1.0+0.81 

GJ 1224 

3000 

-32673 

7 

374 

139 

I.4+O.4I 

GJ 1286 

3000 

-40803 

5 

367 

96 

O.7+O.9I 

LHS 1723 

3000 

42430 

6 

367 

115 

1.0+0.81 

LHS 3799 

3000 

-1641 

6 

602 

74 

I.5+O.5I 


“The error on the absolute RV measurements is estimated to be ~ 100 ms based on 
systematic uncertainties. 

'^The usin^ values reported are the values used in the best fit model. Since values below 
3 km s“^ cannot be confidently measured at our resolution, we adopt a usim upper limit 
value of 3 km s~^ for these stars. 
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Table 6 

Astrometric Results 


Name 


R.A. 

Deck 

Filter 

N ^ 

•1-' sea 

Nf 
' j rm 

Coverage^ 

Years^ 

N_r_bf 

TT (rel) 

TT (corr) 

TT (abs) 


RA. 

V tan 


J2000.0 

J2000.0 







(mas) 

(mas) 

(mas) 

(mas yr~^) 

(deg) 

(km s“^) 

(1) 


(2) 

(3) 


(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

G 99-49 

06 

00 03.52 

-1-02 42 

23.6 

V 

14s 

374 

1999.91-2012.94 

13.04 

4 

191.08+0.97 

2.52+1.57 

193.60+1.85 

307.6+0.3 

97.1+0.10 

7.5 

GJ 300 

08 

12 40.88 

-21 33 

06.8 

V 

14s 

374 

1999.91-2012.95 

13.05 

8 

121.37+0.42 

1.27+0.26 

122.64+0.49 

698.6+0.1 

178.6+0.01 

27.0 

GJ 406 

10 

56 28.91 

-1-07 00 

53.2 

R 

12s 

139 

2000.23-2012.27 

12.04 

6 

413.70+1.61 

1.46+0.17 

415.16+1.62 

4694.5+0.5 

236.2+0.01 

53.6 

GJ 555 

14 

34 16.82 

-12 31 

10.3 

V 

13s 

193 

2000.14-2012.26 

12.12 

6 

161.15+1.45 

0.58+0.22 

161.73+1.47 

682.5+0.4 

331.0+0.07 

20.0 

GJ 628 

16 

30 18.07 

-12 39 

45.4 

V 

10s 

141 

2003.51-2012.58 

9.07 

5 

229.10+2.23 

1.43+0.17 

230.53+2.24 

1191.5+0.9 

185.5+0.07 

24.7 

GJ 729 

18 

49 49.37 

-23 50 

10.4 

V 

11s 

124 

1999.62-2012.75 

13.13 

7 

335.64+1.30 

3.95+0.99 

339.59+1.63 

666.3+0.4 

106.8+0.05 

9.3 

GJ 1002 

00 

06 43.19 

-07 32 

17.0 

R 

9s 

64 

2003.77-2012.87 

9.10 

4 

205.09+3.09 

2.09+0.16 

207.18+3.09 

2034.5+0.9 

203.8+0.05 

46.5 

GJ 1065 

03 

50 44.29 

-06 05 

41.7 

V 

10s 

86 

2003.95-2012.95 

9.00 

5 

99.98+1.91 

1.65+0.24 

101.63+1.93 

1444.8+0.7 

198.9+0.05 

67.4 

GJ 1224 

18 

07 32.85 

-15 57 

47.0 

I 

10s 

170 

2003.52-2012.52 

9.00 

7 

125.04+0.92 

1.50+0.50 

126.54+1.05 

702.3+0.4 

241.0+0.06 

26.1 

GJ 1286 

23 

35 10.47 

-02 23 

20.8 

I 

10s 

135 

2003.52-2012.88 

9.36 

5 

139.19+1.07 

2.28+0.25 

141.47+1.10 

1141.8+0.4 

137.8+0.04 

38.3 

LHS 1723 

05 

01 57.43 

-06 56 

46.5 

V 

14s 

258 

1999.81-2012.75 

12.94 

4 

187.19+0.76 

1.47+0.21 

188.66+0.79 

770.8+0.3 

226.3+0.04 

19.4 

LHS 3799 

22 

23 07.00 

-17 36 

26.1 

V 

9s 

118 

2003.52-2012.70 

9.18 

6 

137.70+1.86 

0.47+0.14 

138.17+1.87 

769.1+0.4 

157.7+0.06 

26.4 


^Number of seasons (Nsea) counts observing semesters where a dataset was taken, and denotes if coverage was ’c’ontinuous (more than one night of data in all seasons) or 
’s’cattered. 

*^Coverage and Years run from the first to last epoch. 



Table 7 

Astrometric Residuals and Errors 


Name 

R.A. Res. Dev. 
(mas) 

R.A. Mean Err. 
(mas) 

Decl. Res. Dev. 
(mas) 

Decl. Mean Err. 
(mas) 

G 99-49 

3.15 

4.70 

2.72 

4.72 

GJ 300 

1.60 

2.58 

2.35 

2.77 

GJ 406 

4.85 

5.18 

5.53 

6.60 

GJ 555 

4.56 

6.18 

4.35 

6.88 

GJ 628 

3.74 

5.30 

4.15 

6.87 

GJ 729 

3.25 

3.43 

4.41 

4.90 

GJ 1002 

0.86 

2.39 

1.54 

3.35 

GJ 1065 

2.01 

3.26 

4.45 

4.23 

GJ 1224 

2.49 

4.60 

2.12 

5.41 

GJ 1286 

1.99 

2.50 

4.28 

4.83 

LHS 1723 

2.38 

3.41 

2.34 

3.97 

LHS 3799 

2.02 

2.67 

2.07 

2.65 


Table 8 

Companion Mass with a 90% Detection Rate 


Radial Velocity Mass Limits Astrometric Mass Limits 


Name 

Stellar Mass 
(Msc/jv) 

3 d 

10 d 

30 d 

100 d 

2 years 

4 years 

6 years 

8 years 

G 99-49 

0.19 

1.0 

1.0 

2.5 

5.5 

23.5 

14.5 

11.0 

9.0 

GJ 300 

0.20 

1.0 

1.0 

2.0 

2.5 

15.5 

9.5 

7.5 

6.0 

GJ 406 

0.09 

0.5 

1.0 

1.0 

2.0 

7.0 

4.5 

3.5 

3.0 

GJ 555 

0.22 

1.0 

1.5 

4.0 

5.0 

41.0 

26.0 

20.0 

16.5 

GJ 628 

0.26 

1.0 

1.5 

2.5 

5.5 

19.0 

11.0 

8.5 

7.5 

GJ 729 

0.17 

1.0 

1.0 

4.0 

4.5 

9.5 

6.0 

4.5 

4.0 

GJ 1002 

0.11 

1.0 

1.0 

3.0 

8.0 

9.5 

6.0 

4.5 

4.0 

GJ 1065 

0.18 

1.0 

3.5 

4.5 

24.0 

31.5 

19.0 

15.5 

12.5 

GJ 1224 

0.14 

1.0 

1.5 

2.5 

5.5 

25.5 

15.5 

12.0 

10.0 

GJ 1286 

0.11 

1.5 

1.0 

3.0 

8.0 

11.0 

6.5 

5.0 

4.0 

LHS 1723 

0.15 

1.0 

1.0 

3.5 

10.0 

13.0 

8.5 

6.5 

5.5 

LHS 3799 

0.14 

1.0 

1.0 

3.0 

3.0 

16.0 

10.0 

8.0 

6.5 
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Fig. 3.— Nightly averaged relative radial velocity measurements are plotted by epoch (see Table H] for a list 
of JD). Vertical dotted lines are used to indicate long time spans between different observing runs (~ 1 to 3 
weeks). After completing periodogram tests on this RV data to seach for periodicity, we find no indication 
of companions around these stars. 
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Fig. 4.— Relative astrometric measurements are plotted over time for both R.A. and Decl. After solving 
for parallactic and proper motion, we see no indications of companions in the astrometric residuals for any 
of these stars. 


31 



































0 20 40 60 80 100 2 3 4 5 6 7 8 

Period (days) Period (years) 

Fig. 5.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for G 99-49. 
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Fig. 6.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for GJ 300. 
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Fig. 7.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for GJ 406. 
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Fig. 8.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for GJ 555. 
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Fig. 9.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for GJ 628. 
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Fig. 10.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for GJ 729. 
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Fig. 11.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for GJ 1002. 
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Fig. 12.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for GJ 1065. 
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Fig. 13.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for GJ 1224. 
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Fig. 14.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for GJ 1286. 
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Fig. 15.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for LHS 1723. 
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Fig. 16.— Fraction of detected companions in a companion mass versus orbital period plot, based on Monte 
Carlo simulations using the RV data (left panel) and astrometric data (right panel) for LHS 3799. 
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